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of integration is the area along the strips. AS an example, 
a rectangular microstrip resonator is analyzed. The accuracy 
of this method is good at high frequencies, but decreases as 
frequency decreases. This method results in small matrix 
sizes (e.g., 4 x 4), shows reasonable accuracy within a 
specified frequency range, and is shown to be flexible since 
it can be used to analyze microstrip lines as well as micro- 


strip resonators. 
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open microstrip line; 3 
shielded microstrip line; 
coupled microstrip line; 


shielded microstrip resonator. 
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(b) 


the case €»=1; 

uhewcascse £1, the electric 

field lines have a break at 

the DAI (dielectric air inter- 

face). 

sketch of the charge distribution 19 
on the microstrip line (taken from 

fez); 

Sketch of the potential distribu- 

tion at the y=h plane (taken from 

[a2 iy. 

LSE (longitudinal section electric) 22 
model for a microstrip line; 

general behavior of the effec- 

tive dielectric constant, €e , 

as a function of frequency, 


obtained from the LSE model. 
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Transmission line models for micro- 33 

strip resonators; 

(a) rectangular resonator; 

(b) transmission line model neglecting 
fringing fields; 

(c) transmission line model with 
fringing fields; 

(d) transmission line model with frin- 
ging fields; due to odd symmetry, 
half of model (c) is sufficient; 

(e) open-ring resonator which can be 


modelled by (d). 
(a) Ring resonator at y=h plane [33]. 36 


(d) magnetic wall model of ring resonator. 


(a) Unshielded microstrip line; Ay 
(b) an equivalent structure 

Equivalent structure for the micro- 46 
strip line used for BEM (Boundary 

Element Method); R is assumed to 

be infinitely large. 

DAI deformed in the vicinity of 4g 
(x;, y,) in order to find the limit 

of dv/ay as (x,y) approaches (X55 h) 

(EQ. 2=11). 

x(¢) discretized to facilitate the 53 


numerical solution of Eqs. 2-22. 


xii 


a 


ao 


-oxD im 10? SRbbm SHEL Nosedimmnert 


gritosizer Lebom enti notesimanest- 
ebisz? anigniat . 


at iw isbom sail moiac imans tt 


DEST Brigniz - 


‘—“izt dtiw isbou ersht fieler saya te 
Cistemays Bho ot of rebisi? ania 


:ineioiTive et (o) fabem to than 


6d oso doidw dozaioset.shitpero — 


»(@Y yd. hefiebom 
Ee) sasig wy ts wetencess gah 


s ToT ahoe st geti7 *o febom [few olfengan 


exis qinseemin bak lebtent 
vutowste jnsfevines as 


(5) - 


ih} 


{e) 
ts) 
(a) 
je) 


(#) 


~otoim edd 262 oxutoyrts dae Lax typ h 
Ciabay ot) WES seh daar. ened aftte 
of bemwane wl A e(beditem treme [¥ 

| on dain ecemittind ‘se 
ke yaiitiniy, od? nt Nenetes Las 
fimit-ad? batt 62 webye a Ps ee 
(at +4) sedoacsaqe (yx) salt ia 


2-6 
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2-10 


Demonstration of ¢-axis divided 


into various segments using 

the data for numerical 

integration. (6 -axis is the path 
defined by the y=h plane in a micro- 


strip line cross-section) 


Potential distribution of a micro- 


strip line on the DAI vs. normalized 


distance with €,=9.9 and w/h=1.0; 
results coincide with those of [12] . 
Surface charge of a microstrip line 
vse normalized distance on the MIC 
with €=9.9 and w/h=1.0; results 
coincide with those of [12]. 
Characteristic impedance of a micro- 
strip line vs. w/h for various €,, 
results coincide with those of [8]. 
Normalized wavelength of a micro- 
strip line vs. w/h for various ; 
results coincide with those of [8]. 
(a) Coupled microstrip lines; 
(b) coupled microstrip lines with 
normalized x-axis. 


Even- and odd-mode characteristic 


impedance of coupled microstrip lines 


vs. w/h for various values of s/h 
and €,=9.0; results coincide with 


[4 Jana [26] . 
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2-12 Even- and odd-mode €, of coupled micro- 
strip lines vs. w/h for various values 


Cee, 


3-4 


of s/h and €,=9.0. 

Even- and odd-mode surface charge 
vse normalized distance on coupled 
mierostrip, ines: €.=9.0, «w/h=1.0 
andiie/N=0. 3% 

Even- and odd-mode potential dis- 
tribution of coupled microstrip lines 
on the DAI vs. normalized distance; 
€-.=9..0.5 (w/h=]..0 sand «s/h=0,.3. 

Even- and odd-mode é&, of coupled 
microstrip lines vs. s/h; w/h=1.0 
and €.=9.0. 

Open microstrip line;, R is assumed 
infinitely large. 

Characteristic impedence of a micro- 
strip line vs. w/h. 

Normalized wavelength of a microstrip 
line vs. w/h. 

Even- and odd-mode characteristic 
impedence vs. w/h for s/h=2.0; 
coupled microstrip lines with 
é€,-=9.0. 

Percent error in final value of é¢ 
vs. Ni Ge is the number of terms 

to which the sum of Eq. 3-24 is 
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w/h=1.5, sAnz0n3} é. =9 0 and 
L,=> (L, is the final matrix size) 


Percent error in Z,vs. Ls (L, is 


' the final matrix size) 


w/h=1.5, s/h=1.0, €.=9.0 and 

Nazto (N,, is the number of terms to 
which the sum of Eq. 3-24 is trun- 
cated). 

A single microstrip line shielded 
on top by a metallic plane. 

Countour for the complex integration 
of Eq. 4-23b; vy, is used to find a 
complimentary kernel. 

&-axis divided into various segments 
for the numerical integration of 
Eqs. 4-23a and 30. 


€&_ of a single microstrip line vs. 


normalized frequency; €,, is obtained 


from quasi-TEM methods. 
ésnoteatsingleimacrostrip line vs. 
normalized frequency; €,.=10.1 and 
w/h=2.53 Ea) is obtained from quasi- 
TEM methods. 

Microstrip line shielded on top by 
a metallic plane; f(x) of Eq. 5-4 
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a= 7 


5-8 


-(x,- d,), and (x; - d,) to 

(x; +d). 

Countour for the integral of Eq. 5=5a. 
MIC divided into several intervals 
for the numerical integration of 

Bq «5-15. 

Rectangular microstrip resonator with 
even mode of resonance. 

Countour for the integration of 

Eq. 5-30. 

Effective dielectric constant of 

a microstrip line vs. normalized 
frequency; &,=11.7 and w/h=0.96. 
Effective dielectric constant vs. 

Nx, (Nx is the number of terms to | 
which Eq. 5-10a is truncated). fh=40.0 
GHzmm, w/h=0.96, &,=11.7 and Mr=4 
(Mr is the final matrix size). 
Calculated effective dielectric 
constant vs. final matrix size, Mr; 
fh=40.0 GHz- mm, w/h=0.96, &,=11.7 
and Nx=500. (Nx is the number of 


terms to which Eq. 5-10a is truncated) 
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B Propagation constant along the microstrip line 


c Dielectric constant of free space 

Ep Mrrectryve retative dielectric constant 

Py Zero frequency e, 

fap Even-mode  «, of coupled lines 

fe 9 Odd-mode™ "ce. of coupled lines 

Ep Relative dielectric constant of the substrate 

(Fig. 1-1) 

lg The guide wavelength of the microstrip line 
& In dielectric 

oP In vacuo 

1 =Ww/(2nh) Pigs: l-Lay lb, and. Ic 

o =(s+w)/(2h), Fig. l-le 

p Charge distribution per unit length on strip 

m Magnetic constant of free space 

w Radian frequency 

a/on Normal derivative (to surface) 


Subscripts and Superscripts 


as Air 

ds Dielectric 

es Even (in case of coupled microstrip lines) 

oO: Odd (in case of coupled microstrip lines) 

t: Transverse component (with respect to the direction 


of propagation) 
tans Tangential component 


nor: Normal component 
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Abbreviations 


BEM: Boundary element method 

DAI : Dielectric air interface 

LSE: Longitudinal section electric 

LSM: Longitudinal section magnetic 

MIC: Metallic portion of the y=h plane in Fig. 1-1 
TE: Transverse electric 

TEM : Transverse electromagnetic 

TLM : Transmission line matrix 

TM: Transverse magnetic 
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I. HISTORICAL REVIEW OF NUMERICAL 
METHODS APPLIED TO MICROSTRIPS 
A. Introduction 
Microstrips are microwave circuits which are constructed 
by printed circuit techniques. Prior to the development of 
printed circuit technology, structures formed of conducting 
surfaces of various configurations were used as electromagnetic 


waveguides. Microstrip circuits have been favoured over 


waveguide technology for many applications, due to their low 
cost and compactness. Also, microstrips are of interest since, 
in principle, the same microstrip circuit can be used for a 
very broad frequency range. Extensive research has been 
carried out on microstrip structures. During the decade 
spanning the years 1945 to 1955, this field of research grew 
so rapidly that a special issue of the IRE Transactions on 
Microwave Theory and Techniques was totally devoted to this 
field [1]. 

The basic microstrip structure consists of a dielectric 
sheet, called the substrate, which is bounded by a ground plane 
on one side and by a thin, metallic circuit on the other side. 
The geometry of this circuit determines the electromagnetic 
properties of the microstrip structure. Significant circuit 
size reduction can be accomplished by using an appropriate 
dielectric material; the relative dielectric constant of sub- 
strates is typically about 10. 

At the heart of many complicated microstrip circuits 
lie three basic microstrip structures, namely, lines, coupled 
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lines and resonators, as shown in Fig. 1-1. These structures 


are analyzed numerically in this thesis. 

Microstrip lines (Fig. 1-1a) are the "second generation" 
of printed transmission lines. The original version, called 
"stripline", was introduced in 1949, [2 ]- A stripline is 
similar to a microstrip line, except that it has an additional 
dielectric and ground plane above the strip conductor. Strip- 
lines are no longer commonly used, since they are harder to 
manufacture than microstrip lines. 

A microstrip line (Fig. l-la) is assumed to extend to 
infinity in the z-direction (i.e., the direction of propaga- 
tion of the electromagnetic waves). In practice, the micro- 
strip line exhibits a characteristic impedance, Z, which is 
used to terminate the line. In this way, the reflection of 
electromagnetic waves is avoided. nus, a Line of finite 
length can be treated in the same way as one that is in- 
finitely long. 

Microstrip lines are normally placed in a grounded 
metallic box having dimensions which are usually considerably 
larger than the line width; the resulting structure is called 
a shielded microstrip line (Fig. 1-1b). By shielding a micro- 
strip line, electromagnetic interference problems can be 
avoided. The effect of these shielding walls on the properties 
of microstrip lines is usually negligible. 

A coupled microstrip line consists of a pair of closely 


spaced parallel metallic strips bonded to the dielectric 
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Open microstrip line; 
shielded microstrip line; 
coupled microstrip line; 


shielded microstrip resonator. 


substrate at the y=h plane (Fig. l-lc). The ground plane lies 
on the opposite side of the substrate. With two conducting 
strips and a ground plane, the structure has two independent 
modes of propagation. 

Microstrip resonators (Fig. l-ld) consist of a metallized 
area in the y=h plane (rectangular, circular, elliptical or 
triangular shapes are commonly used); the dimensions of this 
metallized area are usually of the order of a wavelength at 
the frequency at which the resonator resonates. 

Advances in semiconductor technology have created an 
interest in the integration of microwave circuits, 

Monolithic microwave integrated circuits are active circuits 
operating at microwave frequencies in which the solid state 
devices and the microwave circuit have been integrated. 

In recent years, the terms "stripline", "microstrip" and 
"microwave integrated circuits" have come to be used inter- 
changeably. However, the term "monolithic microwave integrated 
circuits” refers only to “the active circuit described above. 

In this thesis, only microstrip structures are analyzed 

(i.e., microstrip lines, coupled lines and microstrip 
resonators). 

Rapid advances in microstrip technology have stimulated 
much theoretical research; an important goal of this field 
is the computer-aided design of microstrips. The main 
objective of this thesis is to develop computer algorithms 


which may eventually lead to the analysis of a wide variety 
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of microstrip structures. The basic concepts underlying the 
different theoretical methods used to analyze these structures 
are presented. These methods are important for the design 

of microstrips. It is hoped that these theoretical methods 
will lead to the computer-aided design of microstrips. 

In order to analyze microstrip lines, the type of 
electromagnetic waves which propagate in the lines must first 
be determined. A microstrip line resembles a two-conductor 
transmission line and thus, at first, it may appear that the 
structure supports TEM (transverse electromagnetic) waves. 
However, this is true only if the relative dielectric constant 
of the dielectric substrate, €&, tends to unity. Otherwise, 
the line cannot support TEM waves, since the phase velocity 
of the waves in the substrate will be different from that in 
the air. When the operating frequency approaches zero, the 
electromagnetic fields approach a TEM configuration (even when 
E.¥1), which is referred to as the quasi-TEM regime [2]. 

The microstrip can be assumed to be in the quasi-TEM regime 
when the wavelength of the propagating waves is larger than 
the microstrip line width and the substrate height [3]. 

In the quasi-TEM regime, the parameters of a microstrip 
line can be calculated with the simplifying assumption that 
the propagating waves in the microstrip line are TEM. Coupled 
microstrip lines can be treated similarly. Quasi-TEM methods 
are useful not only as a first step towards the exact analysis 
of microstrips, but also as engineering design aids. For 


example, a 50Q microstrip line can be treated as a TEM line - 
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up to approximately 10 GHz. Thus, several quasi-TEM methods 
which are frequently used are explained briefly in this chapter. 
With the continuing advances of semiconductor technology, 
microstrips are being designed for operation at ever higher 
frequencies, where the quasi-TEM approximation may not be 
justified. An exact analysis of microstrips thus becomes 
necessary. Moreover, quasi-TEM approximations can be used 
only for some resonator shapes, such as rectangular resonators. 
In this chapter, some mathematical procedures which have been 
applied to date to the analysis of microwave integrated cir- 
cuits (Fig. 1-1) are discussed. The merits of each method 
(Oe: of 5., the accuracy of the results and the computer space 
and time requirements) are also discussed whenever such data 
have been given in the original papers. However, it is neces- 
sary to emphasize that some of the methods are more analytical 
in nature and that the time spent in formulating them may over- 
shadow considerations of the final matrix size and computer | 
efficiency. Such analytical methods are important if they | 
can be generalized to a wide variety of microstrip structures. 
Thus, a number of methods used by engineers are explained in 
this chapter and the possibility of generalizing them to more 
complicated problems is also discussed. 
Other aspects of microstrips, which deal with higher 
order modes, radiation from sharp edges, microstrip loss, 
surface waves, and box resonances are not included in this work. 


The basic concepts necessary to formulate and analyze 
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microstrips are discussed first. Next, the mathematical 
methods applied by others to microstrips are summarized. 
fhe-methods are classified and discussed according to their 
application to microstrip lines, coupled lines and resonators. 
These methods ore quite diverse and few of them can be 
used to analyze all the microstrip structures shown in 
Fig. 1-l. 

sec. B, titled "Microstrip Lines", is based on the work 
done by Mittra and Itoh [2]. A modified discussion is presented, 
which includes more recent publications. The approach has 
been to avoid mathematical detail, while, at the same time, 
to give a sufficiently comprehensive review of each method 
to communicate its essence. Some quasi-TEM methods are 
presented first. Next, a semi-empirical model is explained 
which describes the frequency behavior of microstrip lines. 
subsequently, the works of other researchers dealing with the 
exact analysis of microstrip lines are presented. 

Sec. C (Coupled Microstrip Lines) is similarly a modifi- 
cation of the work of Weiss [4], together with the review 
of several works published in recent years. 

In Sec. D (Microstrip Resonators), some approximate 
methods for evaluating the resonant frequencies of microstrip 
resonators are explained. Such approximate methods are 
applicable to simple resonator geometries, such as rectangular-, 
circular-, or ring-resonators. Subsequently, other exact 
formulations are reviewed. 


Sec. E (Review Conclusions) compares the methods discussed 
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in this chapter. It is argued that methods involving 
integral equations are better suited for analyzing micro- 
strips than other methods. Furthermore, this section intro- 


duces the methods that are developed in this thesis. 


B. Microstrip Lines 
Quasi-TEM Analysis 
Preliminary Discussion 


For sinusoidal time variations, the electric and 
magnetic field components of the waves propagating in 


microstrip lines must satisfy the steady-state Maxwell's 


equations 
VXE=joul, (1-1a) 
Vx H=fuexE, (1-1b) 
V: E=0, _» thet) 
VH-0, (1-14) 
where 
E is the electric field intensity, 
H is the magnetic field intensity, 
p is the magnetic permeability of free space, 
€ is the dielectric constant of free space, 


is the relative dielectric constant; 
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where f is the frequency. 

In the analysis of microstrip lines, it is assumed 
that the lines are lossless and infinitely long. Con- 
sequently, the fields travelling along the positive z-axis 
are assumed to vary as exp(-jféz). 6 is called the propa- 
gation constant in the z-direction. Eqs. 1, together 
with the assumption that the fields vary as exp(-jfz) 
in the z-direction, result in Helmholtz's equation for 
E and H in the micréstrip line cross section (i.e%; 


the x-y-plane): 


VIEH( ck*-¢*) E=0, (1-2a) 
and 
V H+ (Kk -6 )H=0, (1-2b) 
where 
¥, = (2/9x)°4 (2/ay),, (1-2c) 
ke2mf/c , (1-2d) 


and c is the speed of light in free space. 

Microstrip lines can be fully analyzed by solving 
Eqs. 1-2, subject to the boundary and continuity con- 
ditions for EF and i. However, it has been customary 
to solve Helmholtz's equation for two properly chosen 
potential functions and to obtain the field components 
from these potential functions. This approach will be 
more fully explained in the section "Exact Analysis 
of Microstrip-Lines")(p220). 

In this section, the simpler case of the quasi-TEM 


regime is studied. When €=1, the lowest order solution 
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of Eqs. 1-2 is of the TEM form [2] and [4], with @-k. 
It is then possible to obtain the field components 
from a potential function which satisfies Laplace's 
equation in the cross section 

Viy(x, y) =0. Des, 
where gis the electrical potential distribution in 
the microstrip cross section. 

When E*1, the fields approach a TEM configuration 
for lower frequencies [2d This assumption, as explained 
above, is only valid when f and thus k is small, since 
for €,#1, B-k only as k-0. 

In order to describe the boundary conditions 
imposed on @, two definitions are introduced: DAI 
connotes the dielectric air interface at the y=h 
plane in Fig. 1-1, and MIC refers to the metallic 
portions of the microwave integrated circuits at the 
y=h plane (Fig. 1-1). These acronyms will be used 
throughout the thesis. 

The boundary conditions imposed on #(x,y) are 
as follows (Fig. l-la): 

g(x, 0) =0, (1-3a) 
g(x, h)=v volts on the MIC (1223p) 
The continuity of the displacement vector at the DAI 
implies 
Er BY (¥ W)/ay = ap(x,h*)/ay, ch=82) 


where 
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h=Lim y 
eo a (1-3d) 
y<h 
and 
A=Lim y 


yh (1-3e) 
yh 


Furthermore,—the strip thickness, t, of Fig. l-la 
and 1b is usually assumed to be negligible. With 
this assumption, the free electrons in the strip can 
be regarded as a surface charge layer (x); therefore 
&,09(x, h°) /8y—09( x, h*) /3y=0(x) /s. (1-3f) 
Tneselectricufield, distribution at a. microstrip 
iine.cross Sectionvissillustrated in Fie. 1-2. 
Fig. 1-2a shows the field distribution for the case 
é,=1. . 
When t,#1, due to the boundary condition of 
Eq. 2= 3c, therelectricwiverd lines refract at the 
DAL (Fis. 1=—2b)". 


At this point, the parameters used to characterize 


microstrip lines are introduced fea The capacitance, C, 


per unit length of the line is defined as 
C=0/V, (1-4a) 
where V is the potential on the strip (Eq. 1-3b) and 
Q is the total charge on the strip given by 
Q=f_¥/2o( x) dx. (1-4b) 
When €.=1, the capacitance per unit length is denoted 


by C,.- The characteristic impedance of a microstrip 
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air 


él 


(a) 


dielectric 


Fig. 1-2 Sketch of the electric field configuration at 
a microstrip line cross section 


(a) the case &=1; 


(b)+ethe case €#1; the electric field lines refract 


at the DAI (dielectric air interface). 
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line with €&=1 is defined as 


Zo=1/( eC.) : (1-4c) 
For the case #1, the characteristic impedance is 
Z=2,.( C./C) 1? | (1-44) 


The propagation constant along the z-axis is 
related to the wavenumber, k, as follows 
B=k(C/C,) (1-4e) 
The factor (C/C,) is called the effective dielectric 
constant and is denoted by & . Thus, 
6 6.x, (1-44) 
The physical meaning of € , is as follows: 
Assume that the dielectric substrate of Figs. l-la, 
b, or ec is removed and the whole region around the con- 
ductors is filled with a dielectric material with a 
relative dielectric constant equal to ég, then Z 
and $ will still be given by Eqs. 1-4. However, €, 
is frequency dependent. In the quasi-TEM approximation, 
the value of €é,at zero frequency is obtained. 
Summarizing, in the quasi-TEM method, the Laplace 
equation (Eq. 1-2e) is solved for the potential 
distribution, g, at the microstrip line cross section. 
The boundary conditions imposed on gare described in 
Eqs. 1-3. The characteristic impedance, 2, and the 
effective dielectric constant, Ee, of the line is 
calculated from (Eqs. 1-4). In the following 
sections, the methods used by different authors for 


the quasi-TEM analysis of microstrip lines are summarized. 
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Quasi-TEM Methods 

Wheeler 1 nee conformal mapping for solving 
the Laplace equation at the microstrip line cross 
section. 

The mapping transforms the original microstrip 
line into a parallel-plate structure bounded by two 
vertical magnetic side walls (i.e., a surface on which 
the normal component of Eis zero). The DAI is mapped 
into a curved surface. The curved DAI is then approxi- 
mated by straight lines parallel to the coordinate 
axes. Microstrip lines are classified into three cate- 
gories, according to their width. Different mapping 
functions are used for each category. 

The results obtained by Wheeler [2 | are extensively 
used in microstrip designs, because of their low cost and 
ease of programming. Unfortunately, the method cannot 
be generalized to the exact analysis of microstrips, 
Since conformal mapping is applicable to Laplace's 
equation only. 

In the finite difference (relaxation) method, the 
unknown field distribution is quantized at discrete 
intersections of coordinate grids called the net - or 
mesh-points [5 ] The potential function can be 
expanded in a Taylor series about each mesh-point. 

The potentials at the conductors are known. In 
this way, a matrix equation is obtained which can be 


solved by a matrix inversion. However, when the 
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number of mesh-points is large, the matrix inversion 
becomes very time ene tithe. 
A superior method which reduces the computer time 
for finding the unknowns is the relaxation method [6] 
and [7]. The relaxatia method is a numerical algorithm 
for finding the unknowns by the method of successive 
approximations. 
The finite difference method is easy to formulate 
and program. The method can be easily generalized 
to the exact analysis of microstrips, but the final 
matrix size is usually large (c.f. Finite Difference 
Method applied to the wave analysis described on page 23). 
In the variational method in the Fourier Transform 
domain (s], the strip is considered as an infinitely 
thin layer of charge. Thus, Poisson's equation is 
solved in the cross section of a microstrip line 
(Fig. l-la) 
Vie( x, y) (1/2) p(x) 6(y-h), (1-5) 
Taking the Fourier transform of Eq. 1-5 in the 
x-direction and applying the boundary conditions yields 
(7, b) (7) /fel7|[ 14e,coth( |7|h) ]}. (1-6) 
Eq. 1-6 leads to a variational form for the 
capacitance per unit length, C,.in the Fourier domain. 
The variational method has the property of inclu- 
ding only the charge distribution in the final equation. 


The accuracy of the results depends on the assumed 
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charge distribution. The method does, in general, 
provide a lower bound to the value of the capacitance. 

The typical computer (CPU) time needed is less than 10 
secs. on an IBM 7094 computer [8]. 

The integral equation method uses Green's functions 
and can be used in cases where the strip thickness 
cannot be neglected [9 | and [10] . Define 

G(x, ¥; x. y4) 

such that 

VeG( x, ¥; Xo, Yo) —{ 1/e) 6(x—x,) 6(y-y,). Gee 
G is the potential at (x,y) due to a charge located at 
(x,,» y,)- Applying the superposition principle, the 
potential can be expressed as 

(x, ¥) =f G(x, ¥; x0, Ye) (xe 7.) dy. (1-8) 

The potential distribution is known on the MIC, say 
@sV; Eq. 1-8 is an integral equation with unknown /, 
and can be solved by discretizing the equation. The 
typical computer time needed is 30 secs. for a resulting 
MaLVIxX¥sizer of I3ORX 30. 

The method of moments [12] is closely related to 
the integral equation method. In this method, a 
function similar to the Green's function is found for a 
very narrow but finite line. In this way, the integral 
of Eq. 1-8 is discretized. 

In the generalized Wiener-Hopf method [12] and [13], 
the potential and the charge distributions are found by 


considering their properties in the Fourier domain. 
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This method uses Eq. 1-6 obtained above. 
Taking the potential of the strip as 1V, we can 
express gat the y=h plane as 
g(x, h) =h( x-w/2) u( x-w/2) 
+{ u(x+w/2) -u(x-w/2) J +h( -x-w/2) u(-—x-w/2). (1-9) 
where u(x) is a unit step function and h(x) is the 
unknown potential function on the DAI. Taking the 
Fourier transform of Eq. 1-9 with respect to x we 
obtain 
(7, b) =2[ sin( 7/2) /7] 
texp( jyw/2) H( 7) +exp( -jyw/2) H{ -7). (1-10) 
where H(¥) is the Fourier transform of h(x)u(x). 


Comparing Eqs. 1-6 and 1-10 we obtain 
P( 7Y{el7i[ 1+e,cot h( |7|h) ] } =2sin( yw/2) /y 


texp(jyw/2) H( 7) +exp(-j7w/2)H(-7). = (1-11) 
Eq. 1-11 includes two oan functions, P(?) and 
H(¥), but both can be found by considering the 
properties of the complex functions involved. 
Firstly, -P (Y) has no poles in the entire complex 
%-plane, since /(x) is a function which is nonzero on 
the MIC only. Secondly, H(¥) must have all its 
singularities in the lower 7-plane, since h(x) is non- 
zero for x greater than w/2 only. H(7) and H(-7) are 
thus expanded in terms of their poles which must coin- 
cide with the poles of the expression in the left-hand 


side of Eq. 1-ll. 
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Evaluating Eq. 1-11 at the known zeros of the 
left-hand side posal in a set of equations in terms 
of the unknown coefficients of expansion for H(7). 
These equations are then transformed into a much more 
rapidly converging set of equations. 

sketches of the charge and potential distribution 
obtained by the Wiener-Hopf method [12] are shown in 
Fig. 1-3. The method gives accurate results even for 
very small matrix equations (~5x5), and provides the 
charge and the potential distributions. However, the 
"Wiener-Hopf method is, in general, difficult to formulate 
and cannot be used as a prescription when applied to 
new problems. The generalization of the method to 
coupled microstrip lines and resonators seems very 
dsitieulcr. 

In the mode-matching method [3 |, the cross section 


of a shielded microstrip line is divided into several 


rectangular regions. In each rectangular region the fields 


are expanded in terms of known solutions of the Laplace 


equation with unknown coefficients, The applica- 

tion of the boundary and continuity conditions yield 
the unknown coefficients. The resulting matrix size 

is large (40 x 40 in [3 ]). The accuracy of the results 


oscillates as a function of the resulting matrix size [3]. 
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Fig. 1-3 
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Sketch of the charge distribution on a 
microstrip line (taken from [12] ); 
sketch of the potential distribution at 


the y=h plane (taken from [22] ee 
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Exact Analysis of Microstrip Lines 
Preliminary Discussion 


In the previous section, the quasi-TEM approach 
to microstrip lines and the concept of é,were introduced. 
However, since the wave velocity in the dielectric is 
different from that in free space, it is impossible to 
support a TEM mode in the microstrip line; as the 
frequency is increased, €g@ and consequently Z are changed. 
In the exact analysis of microstrip lines, called 
the wave analysis, the field components are usually 
expressed in terms of a superposition of the TE (trans- 
verse electric) and TM (transverse magnetic) fields, 
which are, in turn, derivable from two scalar potentials 
fis], namely, 
yh) and >) 
which satisfy the Helmholtz equation (i.e., Eqs. 1-2) 
such that | 
E,=j[ («k’-6*) /p) ¥°°)( x, y) exp( -j6z) (Ta1iea) 
H,=i[ (xk*-8?) /8] ™ (x, y) exp( -j62) (1-12b) 
E.=V,.y exp( -j6z) + ou/f) axVyyexp(-jgz) — (1-12c) 
H,={ ace /p) a,x Vy exp(—j6z) +1yexp(-jfz) (1-124) 
where E is the electric field and H is the magnetic 


field; the indices z and t denote the z- and the 
transverse-components. ae is the unit vector along 
the z-axis. The boundary and continuity conditions on 
the field components are obtainable by Eqs. 1-12. The 
lowest-order mode approaches the quasi-TEM mode, as 


explained above. 
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In the following sections, methods used for the 
wave analysis of microstrip lines are summarized. The 
strip thickness, t, is assumed to be negligible, unless 


otherwise specified. 


LSE (Longitudinal Section Electric) Model 


- Based on the field configurations in a microstrip 
line, Getsinger [25] has proposed a model for the 
microstrip line which supports LSE (longitudinal 
section electric) waves, shown in Fig. 1-4a. It is 
easier to analyze this model than the actual microstrip 
line. The solution for the structure of Fig. 1-4 is 
found and yields 

E.-Er-{ Ex—E 00) /[ 14+G( f/f 5) *) (1-13) 
where f€ eo9is the zero frequency limit of Ee 
(obtainable by quasi-TEM methods) and G and f, are para- 


meters that depend on the characteristic impedance. 


Eq. 1-13 shows the general behavior of the dis- 


persion relation (i.e., ‘evs. .f); it is drawn in 
Fig. 1-4b. The slope of &— goes to zero as the fre- 
quency approaches zero. The dispersion curve possesses 


an inflection point at 


f,=! p/v 3G 


and asymptotically approaches ¢r as the frequency 
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Fig. 1-4 (a) LSE (longitudinal section electric) model 
for a microstrip line; 
(b) general behavior of the effective dielectric 
constant, €— , aS 4 function of frequency. 


obtained from the LSE model. 
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goes to infinity. In the following sections, more 


rigorous methods of obtaining the dispersion relation 


are summarized. 


Exact Methods 

The two potential functions 

¥" and yf") 

of Eqs. 1-12 satisfy a two-dimensional wave equation in 
the microstrip line cross section (Fig. 1-1b); thus, they 
can be expanded in terms of known solutions to the Helm- 
holtz equation in dielectric and in air. The imposition 
of the boundary and continuity conditions on E and H 
yields a pair of coupled integral equations [16]. 
These equations may be transformed into a matrix equation 
and solved numerically, yielding the dispersion relation. 
The final matrix size and accuracy depend on which 
boundary and continuity conditions on the six electro- 
magnetic components are chosen and how they are mani- 
pulated. Similar methods are presented in [17] and [18]. 

The finite difference method discussed in the 
quasi-TEM approach can be extended to apply to the 
wave analysis of microstrip lines. 

After properly applying the discretized version 
of the wave equation and the boundary conditions at all 
of the mesh-points, an eigenvalue matrix equation is 


obtained. 
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The algorithm in this method is straightforward, 
but the final eigenvalue matrix is large (e.g., 100 x 
100 or larger [19} ).and the results are less accurate 


than those obtained by integral equation methods fel 


The mode matching method discussed in the quasi- 

TEM approach can also be applied to the wave analysis 
of microstrip lines [20]. The dielectric and air regions 
in the cross section of a microstrip line are divided 
into a few rectangles. In each rectangle, the two 
potential functions are expanded in terms of solutions 
to the Helmholtz equation. The functions are chosen 

so that the boundary conditions on the electric and 
magnetic walls are satisfied. Using the conditions for 
the continuity of the electromagnetic field components, 
the solutions written in the rectangular subregions can 
be matched. In this way, a set of equations for the 
unknown coefficients 18 obtained. The final matrix 
required for reasonable accuracy is 42 x 42,[2]. 

In the finite element method, the differential 
equations and the boundary conditions are transformed 
into a variational expression, which yields the field 
distributions and é&@ when optimized. The microstrip 
line cross section is subdivided into a number of 
arbitrarily shaped triangles; each corner of a triangle 
is called a node. In each triangle, the fields are 
written in terms of known linear functions of x and y 


and unknown nodal values of the potential functions. 
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Since the potential functions are linearized in each 
triangle, it is pansible to integrate the variational 
expression. 

The resulting expression is then optimized with 
respect to the unknown nodal values, yielding an 
eigenvalue matrix equation in terms of é@e. The resul- 
ting matrix size may be as large ae 126 x 126; the 
obtained results are accurate within a few percent [21] . 

In the transmission line matrix method (TLM), the 
microstrip line cross section is subdivided into a mesh 
as in the finite difference method. However, the func- 
tions are not linearized in between the mesh-points; the 
distance between any two adjacent pairs of mesh-points 
which lie along the coordinate axes is considered to be 
a segment of transmission line. The response to an 
impulse applied at a node is calculated by letting the 
impulse scatter in the mesh and reflect back and forth 
a large number of times. 

This method is easily capable of handling losses 
and discontinuities. The dimension of the final matrix 
size is typically 270 X 270. The iteration number is 
200, the required time is 2.16 min. on an ICL 1906A 


computer, and the required memory space is 20k words [22]. 


Open Microstrip Lines 
The effect of the enclosing walls on the properties 


of a shielded microstrip line is usually negligible (by 
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design). However, some mathematical methods are simplified 
when applied to open microstrip lines (Fig. l-la). 

In the integral equation method [23] . the Fourier 
transform of the Helmholtz equation for the potential 
functions of Eqs. 1-12 is taken with respect to x. The 
Fourier transforms of the two potential functions are 
expressed in terms of known solutions of the Helmholtz 
equation. Separate solutions are written in dielectric 
and in air, which satisfy the boundary conditions on 


the ground plane and at infinity. 
The following conditions at the y=h plane (Fig. 1l-la) 


are used to match the solutions obtained for the dielectric 
and the air regions 
E,(x,h”)=E,(x,h”), on the DAI 
BE, (x,h”)=E,(x,h”)=0, on the MIC 
H(x,h)-H)(x,h”)=J,(x), on the MIC 
H,(x,h”)-H,(x,h”)=J, (x), on the MIC 
where E, is the tangential component of the electric field, 
HL and H, are the x- and z-components of the magnetic field, 
and Jy and Ji, are the x- and z-components of the current 


distribution on the MIC, respectively. 


The matching procedure results in a set of equations 
in terms of the current distribution on the MIC. The 
resulting equations are then solved by trying suitable 
functions with unknown coefficients which approximate the 
actual current distributions. This method of expressing 
the unknown fields in terms of the current distribution 


on the MIC is widely used in microstrip analyses. 
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The resulting matrix equation can be solved, 
yielding €4. The size of the final matrix equation 1s 
usually small and depends on the assumed current distri- 
bution [23] . 

In Galerkin'’s method [ 24] the final equations are 
solved in the Fourier transform domain. The tangential 
component of the current distribution along the z-axis, Jy» 
is non-zero only on the MIC (Fig. 1-1a); Ez, is non-zero 
only on the DAI. Thus, the product of J2(x)eand E(x, h) ie 
identically zero. Applying Parseval's theorem, we have 

L235 s(x) ES( x, bh) dxaj/( 2m) [PT (7) B( 7, h) dy =O. Cie) 
A similar equation can be obtained in terms of J. 
Eq. 1-14 can be used as a boundary condition yielding a 
pair of coupled integral equations in terms of ae and 
de - The equations are then solved by expanding apy 
and Je in terms of known basis functions with unknown 
coefficients. The merits of this method are similar to 
the integral equation method described in the previous 


section. 


Quasi-TEM Analysis of Coupled Microstrip Lines 


Preliminary Discussion 
Due to the symmetry of the coupled lines about 


the y-axis (Fig. l-lc) there are two normal modes for 
the structure, namely, the even mode and the odd mode. 
In the quasi-TEM regime, the mode with equal voltages on 


the left and the right strips is called even, and the 
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mode with equal amplitude voltages but of opposite sign 
is called odd. 

since the field configurations are different for 
the two modes, it is expected that there are two values 


for Z and € , namely 


cranes even-mode characteristic impedance, 
Gauet odd-mode characteristic impedance, 
Cee * even-mode é, 


€e,0 ? odd-mode €, 
The numerical methods applied to the coupled 
microstrip lines in the quasi-TEM regime are ‘presented 


in the following sections. 


Green's Function 

In the Fourier integral method, a function 
g(x,y;x"',h) is found, where g(x,y;x';h) refers to the 
electrostatic potential at point (x,y) due to a line 
charge at (xh), pelle In practice, it is easier to 
consider a charge element of nonzero width and finite 
charge density. 

The strips may be regarded as a parallel arrange- 
ment of such charge elements. Using the superposition 
principle, 

9(x, y) =f" e(x, y:x', h) p(x’) dx’, (1-15) 
where f is the charge distribution over the microstrip 


lines. Let R(?7) be the Fourier transform of p(x) and 
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G(x,y;?) be the Fourier transform with respect to x' 
of &(Xsy3x" sh) then Eq. 1-15 can also be written as 

(x.y) =f"G(x, y; 7) R(7 ) ay. (1-16) 
For the normal modes, the potential values on the strips 


are known. The charge distribution can be found from 


Eq. 1-16. 


- G(x,y;%) can be interpreted as the potential dis- 
tribution due to the surface charge 
P(x) + 1/27) exp( -jyx), (1-17) 
since 
G(x. y; 7) (1/27) f_2e(x, y; x’, h) exp(-jyx) dx’. (1-18) 
Thus, G can be obtained. An alternate method for 
finding the Green's function is described below: 

In the bound charge method, the Green's function 
is first approximated by a function g,(x,y;x',h). 
g.(x,y3x',sh) is the same as the actual Green's function 
when the substrate is removed [26] . The error in 
this assumption is a function g,(x,y;x',h) such that 

g(x, y:x',b)=¢,(x, y:x’,h)+g.(x,y;x’,h). (1-19) 
Since Ze is the solution for the structure without 
the dielectric, it will be the solution to the actual 
structure, provided that there is a charge distribution 
at the y=h plane given by 
p.(x: x’) =e(e,—-1) a, Vgo(x, bh; x’, h). (1-20) 
We may place a compensating bound charge 


p(x;x') at y= 
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so that g satisfies the DAI continuity condition. The 

next step is to find a”(x;x"') and go (xsyix'sh). 

Sate given by 

BeCxr yuk ney ater p(x +2) Bx, yrx sncx’, iol) 

where L is the line described by y=h (Fig. l-la). 
Applying the continuity condition at the DAI to 

g yields an integral equation in terms of the unknown f.. 
Having found fk» Eq. 1-21 can then be used to 

obtain Ee: The desired Green's function is eventually 

found from Eq. 1-19 by adding g, to Be° The time 

needed on an IBM 360/67 computer is 0.72 + 21.4Nr, 

where NT is the number of the lines of data produced 

(e.g., 214.72 secs. for 10 lines of data). The computer 

space needed is 64,000 bytes [26]. 


Mode Matching 


This method is the same in principle as the mode 
matching method described under the quasi-TEM analysis 
of microstrip lines on p. 18, 3 , only this time for 


coupled microstrip lines. 


Variational Method 


The Laplace equation can be transformed into the 
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variational expression 

We Sfx (89/8x) *+{ ap/ay) *] dxay, (1-22) 
where S_is.the,. cross. section of the microstrip line. 
W can be interpreted as the stored electric energy. 
Ifa function, ¢ is, found such. that,W. is.minimized,.then 
¢~ satisfies the Laplace equation. 

' The trial field ¢ can be chosen to be the same 
functions as used in the mode matching method above (Zale 
After the necessary boundary and continuity conditions 
are imposed on the solutions, a set of unknown coeffic- 
ients remain. W of Eq. 1-22 is then optimized with 
respect to these coefficients, yielding a set of linear 
equations in terms of the coefficients. The final matrix 
size is about 4x4. The method provides the upper-bound 
values of the capacitance. The results are poor for 
very narrow microstrip lines, but the accuracy increases 
for moderately wide microstrip lines fear This work is 
complementary to that of Yamashita and Mittra [8], nia 
that their results provide lower bound values for the 


capacitance. 


Variational Method in the Fourier Transform Domain 
This is exactly the same method as explained in 
the quasi-TEM analysis of the single microstrip lines 
[28]. For coupled microstrip lines the method should 
be modified slightly to accommodate for the even and 


the odd modes. The charge distribution is identical 
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on the two strips for the even Modes. = or tne odd 
modes the charge distributions on the two strips are 
equal in magnitude but are opposite in polarity. 

The required computer time per data point is 32 
secs. on a Univac 1108 computer. The efficiency of this 
method depends on which function is chosen to represent 


the charge distribution [28]. 


De Microstrip Resonators 


Transmission Line Models 

One of the simplest microstrip resonators is a 
length of microstrip line called the rectangular 
resonator (Fig. 1-5a). Assuming perfect open circuits 
at both ends (Fig. 1-5b), the resonator will resonate 


at frequencies for which 


L=n(A,/2), n=1, 2,°°-. (1-21) 
conanden ne that the guide wavelength is given by 
Age/(trE,), ae 
we obtain 
f=nc/(2L\é,), n=1, 2,---. (1-23) 


Eo is in turn a function of f. The resonant frequencies 
can be found from Eq. 1-21 and the dispersion relation. 
The results obtained by Eq. 1-21 are accurate within 
a few percent, the inaccuracy being primarily caused by 
the fringing fields at the resonator ends. 

A more accurate transmission line model is shown 


in Fig. 1-5c. The end effects are taken into account 


32 


ae A ee ee) 


bio ait rot vesbom 9¥8 itt ish eee ont 

eth sqiute ova ent nS SitAnERILS Sao’ Bit 
wobusdog ag stimoqggo Ste tid, vase 

St @f Jacog SB sag Said sedi bhaaiupss ett ”- oe 
tat Yo yoreiotYte ent ad upon BO.Lt acon papers 
éronemaet of mronorts Bt. mos toss : noicw > abet’ Rolie 
Tee] noiverste Be Te 


: 


~~ 
i 


lk ane 
= : oe oe, a) cee 
2 'IDP 4S LOR Se. bh ee 


= (obo: até no iee ep ceigt ied 7 
stovenckst qiiteoisdn taelaaie eat Te oe | 
es Cyaeternas afF obs bheo- emis diatuoadta 2° : 

ss iveris ambuo sobtisd {asmvagah eladcel iB) OTAONT . .¢ 
atecroest Lliw t*oranesert sit hdeataeitt) ahite avon te 
as thw 40% celndeupen im 


(q#e—t) o* a i= cs ‘Kio ' a * 
7 


ing 
Vd nevis af Asaseletceien at firs ots teat paeren Tea Ny 


(soot) . Cae Dee 


A 


eget) | 8’ Js ‘Cacamnveen 8) > 
dePoaedpott sneer okt .t to Aottonti s _—" at at 

MOPED US tetera th sus brs IHL + ph howe SOR eee 
iti Stemoos | sus IS-1 i pt gd bentst¥e etikaen ott | 
ut Seeush Shr ka aay gested yosavboamt mh .sneotag wat's ) 
) web -tutnngs ty, a ite at 

deotis B . Lebidor Sat =o pienberk aoe A 

atoohte bre od? tok agen 


% in gpa es 
ice ae 


2 


(b) 


nn ae 


(oy (d) 


: COr 


Fig. 1-5 Transmission line models for microstrip resonators; 
(a) rectangular resonator; 
(b) transmission line model neglecting fringing fields; 
(c) transmission line model with fringing fields; 
(d) transmission line model with fringing fields; due 
to odd symmetry, half of model (c) is sufficient; 


(e) open-ring resonator which can be modelled by (d). 
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by adding the capacitors C, at the ends. The values 
of C, are calculated and are available in [29] ana [30]. 

If the rectangular resonator is "bent" so that 
the two ends come close to each other, the open-ring 
resonator shown in Fig. i-5e is obtained [31]. The 
transmission line model for this resonator is shown 
in Fig. 1-5d, [32] « This resonator shows interesting 
neopeteiee [31]. especially for the principal mode 
(ieee, n=1 in Eq. 1-21). For this mode, the resonator 
has an odd symmetry with respect to the z=0 plane; 
Fig. 1-5c¢ can be short circuited in the middle. The 
resulting model is shown in Fig. i-5d. 

The resonant frequency can be obtained by letting 
the impedance at the right and the left at any cross 


section of the transmission line add up to zero 


(Fig. 1-5d). 
21+Z2=0, (1-24) 
where 
Z5=1/( jo), (1=25) 
and 
Z2=Z/t an[ we L/(2C) ] (1-26) 
Substituting Eqs. 1-25 and 26 into 24 we get 
tan[ ové,L/(2C) ] =w0,Z. (1-27) 


Eq. 1-27 yields the odd mode resonant frequencies of 
the resonator. The accuracy of this method is typically 


4% or better [32] ; 
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Transmission-line modelling can also be applied to 
the ring resonator of Fig. 1-6a; the resonator resonates 
at frequencies at which L, the mean length, is an 
integral multiple of Ag. The model gives accurate 
results only when w/L<0.05, [33] - 


Magnetic Wall Models 


The effect of w and curvature in the ring resona- 
tors are not negligible. For this reason, €g as calcu- 
lated for a straight microstrip line is not necessarily 
the same for the ring resonator [33]. Pier ed neCu.Od. 
curvature in the ring resonator is included by regarding 
the resonator as a cavity with electric walls on the 
top and the bottom, and magnetic walls on the sides 
(Fig. 1-6b). The cavity is assumed to be filled with 
the substrate dielectric material. 

Moreover, it is assumed that the variation of 
fields in the y-direction is negligible; the fields 
are TH with respect to the y-axis; Maxwell's equations 
yield 

J, (erkr,) Y,' (ery) J,’ (e,kry) Ya’ (erkr,) 0. (1-28) 
where J, and Y,are Bessel functions [33]: 

This model gives accurate results when w/L40.05, 

[33]. Based on this theory, the mode chart for the 


resonator is obtained in [34]. 
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Fig. 1-6 (a) Ring resonator at y=h plane [33]: 


(b) magnetic wall model of ring resonator. 
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Fringing Effects 


The resonant frequencies calculated by the magnetic 
wall modelling can be inaccurate up to several per 
cent [35] (Cur samUpemtOslLO%—tor the open-ring resonator 
[ 36 | ), since the fringing fields at the resonator 
edges are not negligible. 

' Theseffectwof, the fringing:fields,.can.be taken 
into account by regarding the resonator as a capacitor. 
The theoretical calculation of the capacitance reveals 
an effective change in the resonator dimensions due to 
fringing. Furthermore, an effective dielectric 
constant is found which reflects the nonuniform charge 
distribution in the resonator, which is different from 
the DC charge distribution. Finally, to obtain the 
resonant frequency, a magnetic wall model is considered 
for which the effective dimensions are obtained 
theoretically and is filled with a dielectric material 
with the effective dielectric constant. 

The results of the above theory show good agree- 


ment with experiments [36]. 


Exact Formulation 

In this section, the exact methods which have been 
applied to microstrip resonators are summarized, 
namely, the mode-matching, self-reaction, Galerkin 


and TLM methods. 
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The mode matching method, as described above, can 
be applied’ to the disc resonator. The size of the 
resulting matrix equation is 20x20. With this matrix 
size, the results are accurate to within 1%. The 
computer time on a TELEFUNKEN computer TR440 is 0.666 
secs. [37] - 

_In the self-reaction method as mentioned above, 
the reaction of the field on its own source, should be 
zeroe The ring resonator has been analyzed by this 
method [38]. However, the generalization of this 
method to other geometries seems impossible. 

| Rectangular resonators have been analyzed by 
Galerkin's method [39] which is a generalization of 
the Galerkin method described previously. The typical 
computer time on CDC G-20 computer is 200 secs. per 
structure. More time is required for analyzing narrower 
rectangular resonators. 

in [40], the six component Maxwell equations are 
solved, including the loss in the dielectric and the 
metals by the TLM method. 

This method can be applied to microstrip resona- 
tors of any geometry. The size of the resulting matrix 
is 704 x 704 and the iteration number is 400. The 
computer time on the ICL 1906A computer is 11.26 mins. 


The computer space needed is 46k words. 
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E. Review Conclusions 

In the previous sections, methods for obtaining the 
properties of microstrip lines, coupled lines, and resonators 
were reviewed. Few of these methods can be easily generalized 
so that they will include the three structures of Fig. l-l. 

Since microstrip technology is advancing rapidly, a 
need exists for accurate methods which can be easily pro- 
grammed, which are sufficiently flexible to include all 
three basic microstrip structures, and which are cheap and 
efficient in terms of computer time and memory space. In 
this thesis, novel numerical methods are developed which 
attempt to fulfil these requirements. 

Due to the sharp edges of the MIC, some of the field 
components can be expected to become singular near these 
edges [42] . For this reason, any numerical method which 
does not include these singularities will likely result in 
large matrices. However, since the total energy around the 
edges remains finite, it must be possible to remove the 
singularities. Following this reasoning, it can be deduced 
that methods which employ integral equations are intrinsically 
better suited for analyzing microstrips (since integration 
is performed over the singularity). Many such methods, 
including methods using Green's functions, were discussed 
in the previous sections. In some cases, such as in the 
"Bound Charge Method", it is a major undertaking to obtain 
the Green's functions. 


In order to utilize the interesting properties of 
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Green's functions, while at the same time keeping the effort 
of finding these functions minimal, it is possible to modify 
the integral equations in such a manner that simpler Green's 
functions can be employed. One such method is the BEM 
(boundary element method) [42] and [43]. This method, as 
well as other novel methods, are applied to the analysis of 
microstrips. 

The analysis of microstrip resonators is a complex 
problem. Therefore, simpler structures are analyzed first. 
The methods developed for analyzing single and coupled 
microstrip lines are then applied to the analysis of micro- 
strip resonators. However, each method is useful and applicable 


to engineering design. 


Fe. Thesis Organization 
Chap. 2 deals with the quasi-TEM analysis of microstrip 


lines and coupled lines. The boundary value problem is 
first transformed into an integral equation which is then 
solved by the BEM. Chap. 2 introduces two new concepts. 
First, the thin MIC sheet is included in the integral 
equation and is represented as a charge layer, and the 
singularity in the charge distribution on the MIC is 
explicitly included in the integral equation. Second, the 
integral equation is formulated in such a manner that a 
simple free-space Green's function can be used, regardless 
of the dielectric substrate; therefore, no method of images 


is necessary to find more complicated Green's functions. 
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However, the generalization of the method presents some 
difficulties. For this reason, the integral equation method 
OfeCnap. 2.25 1Urtner modified in Chap. 3. In ‘this 
chapter, Green's functions are generalized and are made more 
flexible so that they can be applied more easily to the wave 
analysis of microstrips. 

In Chap. 4 the generalized Green's functions developed 
in Chap. 3 are employed. Moreover, Chap. 4 introduces a 
new approach to the wave analysis of microstrips. It is 
shown that boundary conditions for Ey, the y-component of 
the electric field (Fig. 1-1), can be totally separated from 
the other components of the electromagnetic fields. Thus, 
instead of using two potential functions Y® ana ph) 
(introduced in this chapter), it is possible to solve the 
boundary value problem for Ey. The parameters of engineering 
interest, such as €, for the microstrip lines, and the resonant 
frequency of resonators, can be obtained from Fy alone. This 
greatly reduces the effort of analyzing microstrips. 

The domain of integration of the integral equation 
derived in Chap. 4 is the DAI. It seems difficult, if not 
impossible, to generalize this method to microstrip resona- 
tors. For this reason, the method of Chap. 4 is modified in 
Chap. 5 so that the domain of integration will be that of 
the MIC. Chap. 5 also gives an example of how this method 
can be applied to microstrip resonators. 


In Chap. 6, the overall conclusions of this thesis are 


summarized and recommendations are made for further research. 
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II. A BOUNDARY ELEMENT METHOD 


A. Introduction 

In this chapter, microstrip lines and coupled lines in 
the quasi-TEM regime are analyzed by the boundary element 
method (BEM). The BEM is a method in which a boundary value 
problem is first converted into a boundary integral equation, 
integrated over the boundaries of a domain. The boundary is 
then divided into a set of elements over which the functions 
under consideration can vary in different ways, as in the 
method of finite elements [42] and [43]. 

The BEM offers important advances over the finite 
element and finite difference methods; the resulting system 
of equations is much smaller, less time is required to 
solve a problem, and the accuracy is greater than in the 
"domain" type methods, i.e., finite elements 42 and 43. 
Moreover, the BEM can be applied to problems with infinite 
domains. 

The classical boundary solution known as the boundary 
integral method is closely related to the BEM [42] , [43] , 
and[44] . 

As noted in Chap. 1, a boundary value problem can be 
transformed into an integral equation by many different 
ways. A good formulation can greatly reduce the effort for 
solving a problem. One possible formulation uses Green's 
functions. However, as Beckenbach [45] has stated, "in 


actual boundary-value problems it is often difficult to 
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determine the Green's function and in many cases it is easier 
to solve the problem directly by numerical methods than by 
the use of the general theory”. 

To analyze microstrip lines in the quasi-TEM regime, 
the Laplace equation for the potential function at the 
microstrip cross section is transformed into an integral 
equation. The domain of integration of the integral equation 
is the y=h plane (Fig. 1-1). In this way, the boundary 
conditions for the thin MIC sheet are imbedded in the 
integral equation and are represented in the form of 
a charge distribution. Moreover, the integral equation is 
formulated in such a way that the free-space Green's function 
can be used, regardless of the dielectric substrate; thus, 
the method of images is not required to find a more con- 
plicated Green's function. Finally, the obtained integral 


equation is applied to coupled microstrip lines. 


B. Formulation 
For a microstrip line in the quasi-TEM regime, it is 
required to solve the Laplace equation 
Vv, 9=0, (2-1) 
in the microstrip cross section (Fig. 2-la). It is assumed 
that the microstrip is unshielded, the ground plane potential 
is zero, and the strip potential is unity. The boundary 
conditions at the y=h plane are as follows: The continuity 
of the potential implies that 
Pla =Ply-n* (2-2a) 


Further, the continuity of the normal component of the 
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Fig. 2-1 (a) Unshielded microstrip line; 


(b) an equivalent structure. 
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displacement vector on the DAI implies that 


&r(99/8y| yn) {39/8 yu*). (2-2b) 
Moreover, on the MIC 
&r(89/89| a7) 1 89/8] you*) =P(x) | (2-2c) 


where p(x)is the charge distribution on the MIC. The micro- 
strip structure in Fig. 2-la is equivalent to that of Fig. 
2D, [1]. The boundary and the continuity conditions on 
the yarn plane in Fig. 2-lb resemble Eqs. 2-2. 

The infinite geometry of Fig. 2-lb is first approximated 
by the finite geometry of Fig. 2-2, enclosed within a circle 
of radius R; the two geometries become identical when R 
tends to infinity. Consider Green's second identity 

Sl PYPy Wr?) dS= 

Scl o( 8¥/an) -{ 8¢/an) ] ac (2-3) 
where S is the area enclosed by the curve C, and 3/an 
designates the normal derivative. In this work gw is the 
potential function, and W is an auxiliary function called 
the fundamental solution. This is similar to a Green's 
function, by means of which @ can be obtained. The procedure 
is explained below. 


Applying Eqs. 2-3 to regions Sy Sos and S3 yields 


Ss,( ep¥ey-Wv 2) ds= 
Sor, [ oC 8¥/8n) -{ 8¢/an) ] dc (2-Ha) 
Ss,( eYv—-Wre) dS= 
Sr, Vodrgic,[ o( 8¥/8n) —( dy/an) dc, (2-4) 
Ss, VY PY") IS= 


Srucgl 9( 8y,/8n) -y( 8¢/An) ] dc. (2-4c) 
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Fig. 2-2 Equivalent structure for the microstrip line 
used for BEM (Boundary Element Method); R 


is assumed to be infinitely large. 
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where \) designates the union of two sets. Thus c,UT, is 


the curve which encloses S,(Fig. 2-2). 


i 
Multiplying Eq. 2-4b by é, and adding it to Eqs. 2-4a 
and 4c results in 
SaUsq( PVPY-WYp) AS+Se,8-( pYPY-Y"y) aS= 
ScUrUrUcgl 9( 8Y/8n) —{ 8y/Sn) } ac 
+Srqlleg\rglle,l 9( 8¥/An) ~Y( By/an) ] dc. ee 
The fundamental solution is assumed to satisfy 
v2v—-5( xx) 3(y-¥9, 3 (2-6) 
regardless of the boundary conditions. Considering Eqs. 


2-1 and 6, the left hand side of Eq. 2-5 reduces to -a,¢; 


where 
Ye ifs (xiv) isin t he vain” 
a= 
e&, if (x,y) is in the dielectric, | (2-7) 
and : 
Pie Xy Ya) . (2-8) 


ieee (xsy; ) is chosen to have a nonzero distance from 
the origin, then as R tends to infinity, the integrands of 
the right hand side of Eq. 2-5 vanish onC,\JC,UcC,UC,.- 
Thus,as R~% Eq. 2-5 yields 
-8 Q=SAT o( 84/8y) VW 2¢/8y) ] y-a* (—dx) 
+S erl p( 09/8) -W{ -O¢/8y) J y-n- dx 
+ fe rl ( 84/8y) —( 89/8y) } yan ( —dx) 
+f Jo -84/8y) WK -99/8Y) | yu* dx, (2=9) 
Considering Eqs. 2-2, Eq. 2-9 results in 
~ aiG= + e,-1) SJ o( 09/8y) Jyent ax—f_W(p/e) [Yl pat dx 
+(&,-1) SJ p( 8Y/Ay) J yn Axt+STE(e/e)[V] yat dx , (2-10) 
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where / is the charge distribution over the upper strip in 
Fig. 2-2. Considering the odd symmetry of y with respect 
to the y=0 plane (Fig. 2-1b), Eq. 2-10 reduces to 


ast &,-1) ff 9] yn'l (89/8y1,-n*) 4 84/8y|y—n*) J dx 


+S Sal p/e) (Hyon*Hy—x*) dx (2-11) 
The solution to Eq. 2-6 is [42 
ya{1/2n) In n= 1/47) In[ (x-x) P+ yy) *], (2-12) 
where 
r=M x-x) #97) ®. (220) 
Thus 
ayp/ay ={ 1/m) { (y-y) /L (xx) PH y-y) “1. (2-14) 


Eq. 2-11 should be modified if (x;,y,;) is chosen to 
lie on y=h, since Eq. 2-14 becomes indefinite as (x,y) 
approaches (x,y;). Therefore, Tz and l, must be deformed 
invo’a semverretewwirth radiue"r> int the vicinity of 
(x3 y; ) as shown in Fig. 2-3; the limit of the first integral 
at the right-of Eq. 2-1). ts.found-as r-tends. to.-zero. 

It can be shown that,if y=h, then 

Sy,9( 89/8n) aCte,S,,9( dY/An) dC=~(é,-1)%, /2 . (2-15) 
Eq. 2-11 should thus be modified to 

(1/2) (e,+1) p= e2-1) JL gly-n*) [ (8Y/8y|,-0*) 4 89/8y|y—n*) J ax 

+S we(P/e) (Vy-n*Hy—n*) dx. (2-16) 

Substituting the results of Eqs. 2-12 and 14 into 16 yields 
+ 1/2) (e,+1) o(x, h”) = | 

(e,-1) f-2 o(x, h*) (b/m) {1/[ (x—x) #+{ 2h) *]} dx 


—f 2/2[ p(x) /(4ne) J nf (x—xy ?/[ (x—xy) #4{ 2h) *] } dx 
. (2-17) 
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Fig. 2-3 DAI deformed in the vicinity of (x;,y,) in 
order to find the limit of ®V¥/y as (x,y) 


approaches (x;,h), (Eq. 2-11). 
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Eq. 2-17 is an integral equation which is equivalent 
to the Laplace equation and the boundary conditions described 
by Eqs. 2-2. To make this integral equation suitable for 
numerical solutions, three additional modifications are 
' required: 

—the fact that p(x) is singular at t+tw/2 must be included, 

=—there must be a changewotrunits to simplify the 
integrals, and 

—the symmetry of the microstrip line about the y-axis 
must be taken into consideration. 

The fact that the energy should be finite in the 
vicinity of the microstrip edges leads to the conclusion 
that p(x) should vary as rat the edges, [ 42] and [46], 
where r is the distance from the MIC edges; thus let the 
charge distribution be explicitly expressed in terms of 
these singularities at x=tw/2;3 

p(x) =ex(x) /M{2x/w)*, (2-18) 
where X(x) is a slowly varying function which describes the 
behavior of the charge fier iaution away from the singular 
edges of the MIC. The factor € is included in the defini- 
tion in order to make the computations easier. An explicit 
representation of the singularity in p greatly reduces the 
numerical effort and leads to a higher computational 
accuracy with smaller matrix sizes [47] . 

The change of units 

¢=2x/w (2-19) 


together with substitution of Eq, 2-18 in 17 results in 
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where 
n=w/(2h), (2-21a) 
C=2x,/W, (2-2) 
©(¢) =p( wi/2, h*), (2-21c) 
$:=8(¢), (2-214) 
and 
x(¢) =x(-we/2). (2-21e) 


Eq. 2-20 can be solved numerically by the usual BEM procedures 
combined with quadrature formulae for integration. 

In the following sections the algorithm for solving 
Eq. 2-20 is presented for special cases and is then generalized 
to include the analysis of single and coupled microstrip 


lines. 


C. Microstrip Lines with §=1 


When ¢,=1, Eq. 2-20 reduces to 


4m1/N=G,+G,, ( 2-22a) 

where 
Gy = —fILx(¢) /M-<?] nl (¢-4) 2H 2/n) "J ae, (2-22) 
Go=/ilx(¢) /M-<?] n( ¢-¢) ®ae (2-22c) 


Note that $;=1 on the MIC. 
Eqs. 2-22 — the same as those obtained by Acton 47 
using Green's functions. The techniques presented by Acton 
47 to treat the singularities are very helpful in solving 
Eqs. 2-22. The integral Eqs. 2-22, when discretized, yield 


a set of linear equations. The domains of integration of 
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Eqs. 2-22 should be transformed to remove the singularities. 
The singularities of Eqs 2-22b occur ac 
Gael, 
The following transformation removes this singularity: 
¢=sin0. 5m (2-23) 
Considering this transformation and the even symmetry of 
the quasi-TEM mode about the y-axis, and using the quadrature 


formulae for integration [4s], we can write 


Gym n/2) BW xan oye) 2H 2/m) I(t yted 2H 2/0) "1} (2-20) 
where 7; 
¢=sin0. 5m, (2-25a) 
and X;, is the value of X at ¢, 
Xy=xXC¢)) - (2-25b) 


Wj and yy are the parameters of the quadrature formulae for 
integration given in [48]. The computer program for 
evaluating G, is called FIRST and appears in App. A. The 
Singularities of Eq. 2-22c lie at 

¢=41, 


and at 
Yu 
To discretize the integral of Eq. 2-22c it is assumed that 
the slowly varying function X(¢ ) is constant at small 
intervals containing ¢;as shown in Fig. 2-4. Thus, equation 
2-22c can be discretized as 
GaP xuffa( 1/MX*) in( C7) Pade. (2-26) 


By the transformation 
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Fig. 2-4 X(8) discretized to facilitate the numerical 


solution of Eqs. 2-22. 
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¢=siné , (2-27) 
the integral in Eq. 2-26 can be rearranged as 
S fai 1/Mi-<?) in( ¢?44) #a¢=f,8#In( sin®9-sin®9,) ®aa 
=Se;7inj [ sin( 6-9) /( 6-9) ][ sin( ( 6+9) /2) /( (848) /2) ] 
[ sin( (1/2) + (6+8) /2)) /( (1/2) + (0+8) /2)) ]} ®ae 
+Je{?ln( 6-9, *da+f,%2In( 6+8,) 2a@ 


+fo{8ln( @+8,-m) #d@- In4) ( @2-9;) (2-28) 
where 
O=sin "ty, (2-29a) 
Og=sin "fy, (2-29b) 
and 
6,=sin*¢,. (2-29c) 


The last step is performed to remove the singularities 


at 
6=8, (0+8,) /2=0, and (6+8) /2=1/2., 
Let | 
A=6-@, (thus e.g. Ap=0,-4,), (2-30a) 
B= 6+8;) /2, (2-30b) 
and 
C= 0+8,-1) /2, (2-30c) 


To use the quadrature formulae for integration, the inte- 
gration domain has to be between -1 and +1, thus the 
following transformation is employed 
| O={ (85483) /2) +¥( ( 82-8;) /2) oars 
Then Eq. 2-28 can be rewritten as 
S§a( 1/8") in( 7) Fag = 
(62-8;) ftinf ( sinA/A) (sinB/B) (sin€/C) ] dv 
+2in|(Ag'®/Ay*1) ( Ba2/B™!) (C,°2/C,) | 
-{ 6+1n4) (42-4;), (2-32) 
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Considering Eq. rs ap the results are presented in 


terms of forms like 


sinA/A and A,4! 


Since as A or A, tend to zero, these forms tend to unity; 


in the programming language, APL, these indeterminate forms 
are defined as unity. In this way, no extra branchings are 
required in the computer programs when Ay) Ans etc., are 
zero in Eq. 2-32. The program to perform the integration 
of Eq. 2-32 is called INT and appears in App. A. 


D. Final Assembly of the Matrix Equation 
The different parts of Eq. 2-22 can be assembled to 


form the final matrix equations. Letting i and j vary over 
all the possible values, Eq: 2-24 results in 
faing [Gi ix}, 
where, from Eq. 2-24 

Gy y= 7/2) Wink [ (eye) BH 2/m) 7 eyed BH 2/) “15. (2-33) 
ix} is a column matrix containing X,, Xgs0-., 
Siler y Has. 2-26 and 31) result in 


[G2] ix}, 


Thus, in matrix form; Eqs. 2-22 can be written as 

iF} =+([ G1} + [G2] ) fx}, (2-34) 
where {F} is a column matrix, all elements of which are 
equalto 4i/m (Eq. 2-22a). The program FNL in App. A forms and 


solves the matrix Eq. 2-34. 
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The capacitance per unit length of a microstrip line 
with €&=l is the total charge per unit length divided by 
the strip potential, which is assumed to be unity. To find C 


the total charge, p must be integrated over the strip, thus 
C =f, pax=enfix/M ae, (2-35a) 


or 


Com en) £ Wy xy (2-35b) 


in the) computations, €yi1s found inwunitstof €¢ - 


E. Microstrip Line with €,41 
In the case where &#1, Eq. 2-20 must be solved. Let 


the first integral on the right hand side of Eq. 2-20 be 
designated by Ge 
: G5= S28(¢) /[ ($4) #4 2/n) J ae ( 2-36a) 
@ is known on the MIC and is equal to one volt. Thus, 
the portion of Ge which is integrated over the MIC, Ge, 
can be easily found as 
G_=(n/2) tan“ (1-¢,)n/2]{n/2)tan“[(14¢,)n/2] (2-360) 
Using the quadrature formulae of integration [4s], the 
portion of Gy which is integrated over the DAI, Gos can be 
approximated by 
Ga=P SM U/L (Sy) PH 2/n) 7] 41/1 (ttt) 2(2/n) 21, ( 2-36c) 
where 
Cyt , ( 2-364) 
j 


W; and % are tabulated in [48]. G, is the sum of G 


and Go 


G; =G,+G_. ( 2-36e) 
The results obtained for the case ¢€,=l1 can be used for 


the case &#1 when slightly modified. Considering Fig. pans 
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Fig. 2-5 Demonstration of ¢-axis divided into various 
segments using the data for numerical inte- 
gration. (€-axis is the path defined by the 


y=h plane in a microstrip line cross-section). 
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G, can be greater than unity for the case £41; G; can 
lie on the DAI, since the integrals on the DAI do not vanish 
(Eq. 2-20) for #1. The matrices [G,] and [G,] should be 
modified to [3] and [c,] respectively, when Cue 


ThecvLlinialemacrix equation is 


[G,] +1 Ge), ix} 


[ Gg] +f GJ) 5 6} ( 2-37) 


where 

if@fis a column matrix containing the unknown g's 
on the DAI (Fig. 2-5). 

[e.] and [c, | are the same as defined above, (Eqs. 2-33 

to 34) 

(ce; | iSrayModirication of [c,] when ¢.7 1, 

[c, | is a modification of [G,] when ¢,>1. 

[G,] is defined by Eqs. 2-36, and 

{Fy is a column matrix, the i' th element of which is 
~2n( e,tl) % 

Eq. 2-37 should be further rearranged since ?, is known 
only when ¢, lies on the MIC. ‘FY contains some ynknown ?; 
and {$} contains some known G3 rearranging Eq. 2-37 results 
in 

ix} 
{F} =f G]4—-- 
i$} (2-38) 
where {F} is a column matrix all elements of which are 
equal to -2%(€ +1). By solving Eq. 2-38, the charge dis- 


tribution on the MIC, the potential distribution on the DAI, 
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and the capacitance per unit length, C, can be calculated. 
Having {x},the charge distribution can be obtained from 
Eq.e 2-18 and the capacitance per unit length, C, can be 
calculated from Eq. 2-35a. Moreover, {$} provides the 
potential distribution on the DAI. 

Having C, and C, it is easy to determine &, and Z from 
Eqse 1-4. Calculated results are presented in Figs. 2-6 
and 7 for the case &=9.9 and w/h=1.0, and for several 
values of €, in Figs. 2-8 and 9. The results coincide 
with the results obtained by Yamashita and Mittra [8] » and 
Mittra and Itoh 12. The overall matrix size is 15 x 15. 
The CPU time per run on the U. of A. AMDAHL 47/V6 computer 
is 0.08 sec. The computer programs, together with an 


illustrative example are presented in App. A. 
As Figs. 2-6 through 9 show, the obtained results are 


identical to those obtained in [8 | and [12]. However, the 
present method is easier to formulate, is more economical, 
and provides more information than other methods. 
Specifically, the variational method [8 ] provides €e and Z 
only; the computer time is 10 secs. on an IBM 7094, 

The integral equation method fen provides P(x), Cp.) end 

Z and the typical computer time is 30 secs. with a final 
matrix size of 30x30. The Wiener-Hopf method [12] provides 
P(x), O(x,h), €, and Z; the final matrix size is less than 
10x10, but it is essentially an analytical method which 


seems impossible to generalize to coupled lines. 


F. Coupled Microstrip Lines 
Eq. 2-17 can easily be generalized to the case of 
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Potential $¥(x,h), Volts 


Fig. 2-6 Potential distribution of a microstrip line 
on the DAI vs. normalized distance with ¢,=9.9 


and w/h=1.0; results coincide with those of [12]. 
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Fig. 2-7 Surface charge of a microstrip line vs. 
normalized distance on the MIC with €= 9.9 


and w/h=1.0; results coincide with those of. 


fiz2): 
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Fig. 2-8 Characteristic impedance of a microstrip line 


vs. W/h for various €,; results coincide 


with those of [a]. 
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Normalized wavelength of a microstrip line 


vs. w/h for various €,; results coincide 


with those of [8]. 
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coupled microstrip lines. The integration domain of the 


second integral on the right hand side of Eq. 2-17 should 
be modified to include both strips. The formulation 
details are given only for the case €,=l since the procedure 
for €#1 closely follows that of the single microstrip 
line. 

Pop .coupledsmicrostrip, Lines with ¢;=1, Eq. 2-17 


reduces to 


4nbef a, [e(x) /e] inf (x—x) °/1 (xx) *H 2h) “T$ dx. (530) 
Let 
Pea (2-40a) 
define (Fig. 2-10) 
a={s+w) /(2h), ane 
and 
n=w/( 2h) , .  (2-40c) 


Then Eq. 2-39 can be rewritten as 


amb f4ERB ppl €) /e] inf (€-€,) */1 (€-) *+(2/n) "14 ae 
+o" pal -€) /e] Inf (€-€) °/T (€-€d PH 2/n) “13 8, (2-'1) 

where #, and # are the charge distributions on the MIC on 
the positive and on the negative €-axis. For the even 
modes 

Pp §) =pal -€); (2-42a) 
For the odd modes 

pp( £) pal -€). (2-420) 
Considering the singularity of the charge at the MIC dees 
f,(&) can be written as 

py=ex( o-€) /M-{ (o-€) /n]?, (2-43) 


where X is a slowly varying function of ¢. 
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(b) 


Fig. 2-10 (a) Coupled microstrip lines; 
(b) coupled microstrip lines with 


normalized x-axis. 
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For the even modes it is assumed that the voltages on 
both strips are +1 volt and that the voltage is zero on the 
-ground plane. For the odd modes, it is assumed that the 
voltages are +1 volt on the right strip, -l1 on the left 
strip, and zero on the ground plane. 

The point és is Chosen to lie on the right strip; thus 

é,=1. (2-44) 


Combining Eqs. 2-41, 43, and 44 results in 


an=szint x (o-€) /Ma-€) /nl 4 
{ant ( £8) °/((€-) 24 2/m) 9] sta (E+E) */( (E48) H2/) DEAE 


where the plus sign holds for the even modes and the minus 
Sign for the odd modes. 
The transformation 
¢=(0-€) /n (2-46) 
simplifies Eq. 2-45 to 
ann SA x(t) /Mi<*) Lal (41-4) ZL SX *H 2/0) 71] 


#in[ ( 1’ %)®/(( ty Y)*#H(2/n)) 13.45, (2-47) 
where 
os =2(0/n) Xs (2-48) 
Eqe 2-47 can be discretized as follows 
where 
Gy yl¥e in tt) 7/M Fae, (2-50a) 
Ga y=/igi in( ty) */M Fae. (2-500) 
Gs y=-Winl (6-4) 74 2/n)*), (2-50c) 
Gy yon (6) 4) 24+ 2/n) 7). (2-504) 


X; was defined in Eq. 2-25b, Wy and 4; are the parameters 
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of formulae for quadrature integration available in |48]. 

Eq. 2-49 can be written in matrix form as 

{ U} =([ Gy] +f Gel) ix} +([ Gs] +f G,] ) Ex}. (2-51) 

where {U} is a column matrix, all elements of which are 
equal to 41, and {X}is a column matrix with elements ae , 
The plus and minus signs refer to the even and the odd 
modes, respectively. To evaluate Gyajy» tet 
f=sind, (2-52) 


Then, it can be shown that 


Gy; =So3"In{ [ sin( ( 8-9.) /2) /( ( 8-8) /2) ]® 
[ meos( (8+8,) /2) ]°/[ ( (1/2) + (048) /2)) ( (1/2) + (6+8,) /2)) 13} a0 
—2( 82-8 ;) In( 27m) +2 f,12in( @-9,) dé 
+2 fg{in[ 1 6+9,) ] dO +2 f92in( 1+9+48,) da. 
(2-53) 


Goi; can be found Similarly. 


Coupled Lines with €,#1 


For coupled lines with Conf legs o2—1 7 should be 


modified as 
—2n( +1) 8=2( 8,-1) f28( €) /T ( &-€s) +4] dé 
+f x(0-€) /M{ (a-€) /n) 7] {In €-€) 7/L ( €-€) 2H 2/n) 7] 
+In( §+§,)7/[ (€+€) *4(2/n) “]} dé, (2-54) 
&( £) =6( né) =8(x). 
The + and - signs in the right hand side of Eq. 2-54 


where 


refer to the even and odd modes, respectively. fo 


simplify the integrations we introduce the following 


change of variables 


C=(G-E)/y , (2-55a) 
O= 6/y ] (2-55b) 
¢. =26—-¢; . (2-55c) 
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Consider the integrals of Eq. 2-54. The potential is 
unknown on the DAI and the charge is unknown on the MIC. 
Thus, the integration domain is divided into three 
regions (i.e., the DAI in between the two strips, the 
strips and the DAI outside the strips) to make the 
numerical integration possible: 
o' -1)v/2+0' +1) /2; 1<¢<o’ 
¢=4sinO. 5nv; =1<¢<1 
+147); —oce<—f (2-55d) 
Substituting the change of variables for Eqs. 
2-55 into 54 and following the same procedures as in 


the previous sections leads to 
{F} {[E.] [Ea] ) £8} 4 [ Es) +{ £4] ) 1} 

+ [Es] +[ Ee] ) x} +([ Ex] +1 Be] ) ix$ (2-56) 
where {#} is a column matrix which contains the unknown 
values of the potential on the DAI in between the 
strips. Similarly, {3} is a column matrix with the 
unknown potential values on the DAI outside the strips. 
{X} is a column matrix representing the unknown charge 


distribution on the strips. 1FY is a column matrix with 
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elements -21(€,+1)$;, if. the point i lies on the DAI and 


-21 (+1) OP 2060-1) Sire VeH Be Prat dg , 
if the point i lies on the MIC (¢.f. the first integral 
in the right hand side of Eq. 2-54). 
LE} ] through [Eg] are defined as 
Ey, ytEe, yt (2-1) (o" —1) /n] Wy 
11/0 (Sat 24( 2/n) 7) 41/0 (6 Sp) 24 2/n) 713; 
—wS¢,<g'; 1<¢,<0'. (2-572) 


at - econeag ant ies 

OEM on? no Neer et ote ws 
esutt oii pai eis ‘ot a si hs Bi 

wit weqeare ow SAT neonate 3 oat ua ‘ea waa 7“ 

a7 agpm oF Ce ite one watetiro bas ait ‘eis ety 


| ws eeldeliny Te egeiids eis ap.rmnni dui e 

ABR De = aiian asist sehiwod Tok ite, #2.citead ie 

oF ae ol ehoifies 2a fred er _ 

ee ee ay 

(ae - ix) (Le Se eC Gee Sa 
qvedunuy edt eciateas dataw x bata ei LOD: f ae £85 cnt a 
ait caenied af TAT sag mo facinetag, oat he esitey . 


iv ew xz eos um pote Pos a ae jay ec velt Siebel £ : 


>On iy ne ria: ees 
a4, 25 o Ty bey 16 fe rie DAT Ody 14 Bers 


sais Hwowsirur syer yobtnees TPP sees 
2 . a * re ; 


‘ 
he ai ua 


ws. tA. eur -n6 aga t ha 
oe bas ie ee iby 
leqastid teat ede 3S 2h Uae Lr Be whe fave id 
a “es 4, ay 7 At atte ve om 
| be * shew Va mg 


al 
rh 


a5 “pead. 38h & 


s UP f 
ee as . Wy et J : : 


a ey 
i. au 1 


/ \a~ , : : : De, 
. ot fe ap Cm ie 


. 69 
Eg ytE, yal 2(e,-1) /n] Wy 


1/7 (Ss-¥p BH 2/n) 7) 41/0 (6s Sp 27H 2/n)*]}; 
SE Sg! 5 S$ ,<—1; (2-57b) 
Es, ytEg, y=ml,/2- 
{ In[ (¢s-¢)) 2+( 2/n) *] snl ( 6-¥) 24+ 2/n) 7]; 


—oSfy<q' 5 S$ <1 (2~57¢) 
Ey, ytEs, gE 1/M-¢?). 
[In( ¢s-<) ?+In( ¢;' -¢) 7] dé 
| SE So"; 156,51 (2-574) 


Moreover, different numerical procedures are 
required in evaluating [E. | and [ Eg] depending on which 


region’ the point 6elies tin > *" Let 
Fay if 1<¢,S0" 


Ee, 4 Gy y if —1 <¢,51 


(2-58a) 
Ay y if —o<¢,<—1, 
Thus 
[Fi] 
[Ee] =|L6) 
[H,] (2-58b) 


Finally, considering Eq. 2-56, some of the $;'s 
appearing in iF} are unknown and should be transferred 
to the right hand side of the equation. Eq. 2-56 
can thus be rearranged as is 
(st 

(Feb -TGeHts) (2-59) 
for the even modes. Similarly a matrix equation can 
be obtained for the odd modes. 

By solving the final matrix equations the even- 
‘and odd-mode characteristic impedance and ¢, can be 


evaluated. The computer programs are presented in 


App. B. The results obtained by this method are 
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presented in Figs. 2-11 and 12, and coincide with the 
results in [4] and [26]. The charge and the potential 
distributions at the DAI are calculated and presented 
in. Bigs. 2-23 and 1+. 

Fig. 2-15 shows the €, vs. s/h obtained for the 
case 


w/h=1.0, and €éy; =9.0 


The final matrix size is 35x35 and the CPU time on an 

AMDAHL 47/V6 per point per mode is 2.165 secs. For the 
sake of comparison, we note that with Green's functions 
126], the required time on an IBM 360/67 is 22.12 secs., 
to obtain the same results hein the same accuracy; the 
variational method ]28] provides only €, and Z with a 

required time of 32 secs. on a Univac 1108, and is less 


accurate (within 2%) than [26]. 


G. Summary 


In this chapter, microstrip lines and coupled lines in 
the quasi-TEM regime were analyzed by the BEM. Laplace's 
equation at the cross section of a microstrip line 
(Fig. l-la) should be solved for the potential dis- 
tribution 

V*y=0, (2-1) 
The boundary and the continuity conditions at the y=h plane 
were 
P| year =Pl yau* (2-24) 
On the DAI, 
&,( 89/8Y|yn-) {99/8 yan*) (2-2b) 


and on the MIC, 
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Figs i2a1) 


0.0 1.0 w/h 2.0 


Even- and odd-mode characteristic impedance 
of a coupled microstrip line vs. w/h for 
various values of s/h and €,=9.0; results 


coincide with [4 | and [26]. 
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Fig. 2-12 Even- and odd-mode €, of coupled microstrip 
lines vs. w/h for various values of s/h and 
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Fig. 2-13 Even- and odd-mode surface charge vs. 
normalized distance on coupled microstrip 
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€,(89/8F\yen-) T 89/8y| pant) = P(X) /e (2-2c) 


where (x) was the charge distribution on the MIC. 
Green's second identity was used 


Ss( ev? 29) dS= 
Sel o( 8¥/an) —( ag/an)] ac (2-3) 
where y was the free space Green's function. Using Eq. 2-3, 
both Laplace's equation and the boundary conditions were 
transformed into an integral equation on the y=h plane. 
+ 1/2) (e,+1) o(xy h”) = 
(@,-1) f= o(x, bY) (b/n) {1/[ (x-xy *H{ 2h) “] } dx 
~f 2/3 p(x) /(4me) Jin{ (xx) ?/[ (xx) #2) *4]}4x (917) 

The singularity in the charge distribution was 

inserted explicitly into Eq. 2-17 

p(x) =©x(x) /M{ 2x/w)?, (2-18) 
where X(x) was a slowly varying function. This reduced 
the final matrix size and increased the accuracy [59]. 

The simple case of a microstrip with ¢,=1 was analyzed 
first. The singularities in the resulting equations were 
removed using appropriate transformations. The case &,#1 
was then analyzed. The numerical results obtained for the 
single lines were shown in Figs. 2-6, 7, 8, and 9 and 
coincided with results obtained in [8] and [12]. 

The method was then applied to coupled microstrip 
lines. The results obtained by this method were presented 
in Figs. 2-11 and 12 and coincided with the results obtained 
in [4] and [26]. Additional original data were presented 
in Figs. 2-13, 14 and 15. 
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Conclusions 

The BEM has been successfuly applied to microstrip 
lines and coupled lines in the quasi-TEM regime; the analysis 
is capable of evaluating the characteristic impedance, Ez. , 
the charge, and the potential distributions. The only other 
method yielding all of the above parameters and distributions 
is, to the author's knowledge, the Wiener-Hopf method [12]; 
however , it seems impossible to generalize the latter method 
to coupled lines and resonators. 

The singularity in the charge distribution is included 
in the formulation. Moreover, the BEM has been modified by 
using the quadrature formulae for integration. This has 
resulted in increased accuracy and in a reduction of the 
size of the final set of equations. For microstrip lines, 
the final matrix size is 15x15 and for coupled lines, 
wher firialematri xsi zes 1s135x95. 

The weaned is very efficient and accurate; the CPU 
time per run on the AMDAHL 74/V6 computer is 0.08 secs. for 
the analysis of microstrip lines and 2.165 secs. per run 
per mode for the analysis of the even- and odd-modes in 


coupled lines. 
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III. AN ALTERNATIVE INTEGRAL EQUATION METHOD 


A. Introduction 

In Chap. 2 the BEM was applied to microstrip lines in. 
the quasi-TEM regime. It is possible to generalize this 
method to the wave analysis of microstrip lines, but the 
fundamental solution will be difficult to obtain. 

To solve the Laplace equation, Green's functions are 
customarily defined so that they satisfy Poisson's equation; 
the forcing function is chosen to be a point charge. In 
this chapter, instead of Green's functions, a fundamental . 
solution, ~, is chosen so that it satisfies the Laplace 
equation. Moreover, the boundary conditions on YW are 
-defined such that the final integral equation is expressed 
in terms of the charge distribution alone. The method leads 
to expressions for the characteristic impedance and the 
effective dielectric constant for microstrip lines. The 
method is then applied to coupled microstrip lines, and in 
the following chapters it is generalized to the wave analysis 


of microstEuips. 


B. Formulation 

For a microstrip line in the quasi-TEM regime, Laplace's 
equation for the potential function y should be solved in 
the cross section of a microstrip line. The boundary con- 
ditions are given in Eq. 2-2 and are demonstrated in Fig. 
3-1. The radius of the semicircle, R, in Fig. 3-1 is 


assumed to be infinite. 
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Fig. 3-1 Open microstrip line; R is assumed infinitely 


large. 


(2 


7 


i) 


it oe 


F A 


vistiaitne béwmide 61 A oni l qitteowotm New 


a: 


wee 


fy Ly 
aa " 


” 

ae A . - 
‘ft £ : att Tah 
ot ei - 


a 


ao 


| 
¢ 


Green's second identity can be applied to S, and 8, 


of Fig. 3-1. 


Ss,( py Wi) ds= 
ScUregWe,[ p( 64/8n) —-{ 4¢/an) } dc; (3-1a) 

Ss, py? ¢) dS= 
Scbrel ( 84/8n) ( 8¢/an)] dc. (3-1b) 


where, as shown in Fig. 3-1, the region Sy is enclosed 
by paths Ci» rj : Ce, and C,; the region S, is 

enclosed by paths le and Co3 g is the desired potential 
distribution; ¥ is an auxiliary function (fundamental 


solution), and 9/gn designates the normal derivative. 
Adding Eq. 3-la to b results in 


S3Vs( 9Vy—-WV7¢) ds= 
Ser Pesle,[ o( 8¥/an) {3p /An) ] dC+feur[ o( 8Y/an) —( 8y/én) ] AC. 


(3-2) 
~ is zero on ( (Fig. 3-1). Wis also defined to be zero 


on C,. Moreover, g and a¥/an tend to zero at infinity. 
Thus by letting R tend to eRe Eq. 3-2 reduces to 
Savs,( pVY-W7¢) dS= 
Srl 9( 8¥/8n) —( d9/8n) ] dC+f;,[ o( 8¥/8n) -{ Ay/an)]dc, (3-3) 
where it is required that 
Vie=0. (3-4) 
Up to this point, the procedure is almost the same as 
that in Chap. 2... However, the auxiliary function W is 
assumed to satisfy Laplace's equation 
V7 ~=0. =>) 
This choice makes the final equation much simpler and 
results in a very simple numerical procedure. Moreover, as 


will be explained in the next chapter, it is easier to 
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81 
generalize this method, rather than the one in Chap. 2 to 


bine wave analysis of, microstrip lines. 
Considering Eqs. 3-4 and 5, Eq. 3-3 can be written 
as 
O=ff 9 aY/8ylya7) V8 9/A8y|y47) J (-dx) 
+f J 9( -8Y/8yly-0*) 99/8 yly-n*) | dx 
=f J 9 6Y/ay|y-a7) -9( 8Y/8y| yn") | dx 
fT 89/aylya) -H( 8¢/8yly2*) ]dx © (ane) 
By properly imposing the boundary conditions on W, Eq. 3-6 
can’ be reduced to a simple integral equation in terms of 
the charge distribution alone; let 
(Aya) =er(Hyatds (3-7a) 
and 
( OY/Ay| ya) 1 8Y/Ay|y-n*) =O. 5[ 6(x-x) +5(x+x,)]. (3-7b) 
Substituting from Eqs. 3-7 into 6, and considering the 
boundary conditions on @ at y=h (Eqs. 2-2) the following 
Take ed equation is obtained 
O=p( xs, h) —Sme(P/&) (Yly-n*) dx. (3-8) 
To find the function wy which satisfies Eqs. 3-7 let 
(considering an even symmetry in a microstrip line), 
JoA( 7) cos( yx) sinh(yy) dy; 0<y<h 
y= 
(7) cos( yx) exp(—-yy) dy. h<y. (3-9a) 
Was defined in Eq. 3-9a satisfies Laplace's equation. Imposing 
the boundary conditions of Eqs. 3-7 on wpyields 
A(y) sinhyh=e,B() exp(—yh), (3-9) 


and | 
B(y) {2/m) exp( yh) cos(yx,) /{7[ite,/tanh(yh)]}.  (3-9c) 
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Substituting w from Eqs. 3-9 into 8 results in 
SSurgt P=) 72] cos(yx) cos(yx,) /{y[ 1+e,/tanh(yh) ]} dxdy=n/2. (3-10) 


where x;is chosen to lie-on the MIC. 
Eq. 3-10 is solved for the single microstrip line and 


coupled microstrip lines in the following sections. 


C. Single Microstrip Line 


As in Chap. 2, the singularity in the charge distribu- 
tion is explicitly included in Eq. 3-10 for reducing the 
final matrix sizes 

p(x) =ex(x) /M~{2x/m)*, (3-11) 
Eqe 2-10 can be written as 


Sot [Ae x(x) cos( yx) /M-( 2x/w) *dx} 


cos(yx,) /{7[ 1+e,/tanh( yh) ] } dy=7/2. (3-12) 
To solve Eq. 3-12, let 
T=f_/2x(x) cos(7x) /MA 2x/w) 24x: (3-13) 


"oneiecine the even symmetry of the quasi-TEM mode ina 
microstrip line, the slowly varying function X( x) can be 
expanded in a Taylor series as 

X(=) Rotkixte +, | (3-1!) 
Thus 


T=f_W/R( K54K,x*+ «+ )cos(yx) / Mi 2x/w) *ax. (3-15) 


Applying the transformation 


sin @=2x/w 


-to Eq. 3-15 and using the tabulated integrals of [ 49] yields 


1-4 0/2) [KTo( 0. 57) 4KyI o(0. S70) +--+]. (3-16) 


Substituting from Eq. 3-16 into 12 yields 


£ KS5I (0. Syn) cos( 7x) /{7[ 14¢,/tenh( yh) ]}ay=1. (7-2?) 
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The following transformations are further used to simplify 


Eq. 3-17 


a=yh, (3-18a) 
t2x,/w, 7 ( 3-18) 

and 
n=w/(2h). (3-18¢) 


Applying Eqs. 3-18 to 17 results in 

E KaSel ex( na) cos( ant) /[a(1te,/tanha)]da=t, (3-19) 
Choosing N points for ¢€, and truncating the sum to N-1 
results in an NxN matrix equation from which the charge 
distribution on the MIC can be found. However, for most 


practical values of y the only «term of,significant value is 


Kote this 
KS6I o( 12) cos( ant,) /a( 142,/tanha) ]dawt. (3-20) 
Onlysone).g;.1s;needed,insorder.toj obtain, Kein Eq. 3-20, 
so let 
{,=0, (3220) 
and hence 
K.Sod o( 7@) /[ a( 1+e,/t anha) J dasi. (3222) 


A direct numerical integration of Eq. 3-22 leads to 
inaccurate results due to the oscillatory nature of J,(n4); 
neither is the integral listed in available handbooks of 
integrals or Hankel transforms [49] and [so]. Thus, the 
fraction in the integrand of Eq. 3-22 is expanded as 
- 1/(1te,/tanha)={1/( 2,41) ] (1-04) E (-1) 5e eae ( 3-23a) 
where 


r+ e,-1) /(e,41). (3-23b) 
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Substituting from Eq. 2-23a into 22 yields 
KE, (—1) "Sal (na) [ee 4) /q daske,t). (3-24) 
Using the tabulated integrals of [49] and manipulating, 
Eq. 3-24 results in 
KZ, (7) "inf [2(n+1) +\4(n41) 247,°] /( 2nt4n*in®) }(e,41). (3-25) 
The capacitance per unit length can be readily found from 
BGs 3-25. For. €:=1 
| C.-2/{ eln[ (24+44n%) /n]}, (3-26) 
andstor { é€.41 
Co{ e,+1) /{ ef, (1) Sn[ [2( att) +4) 4497] /( enn) Jf. 
(3-27) 
Considering Eqs. 3-26 and 27, €e is given by 
£o-f (e,+1) /2] Inf (24-447”) /n] / 
{3 (—7) inl [ 2( a1) +4(n41) 2497] /(2nt4n24n*) ]}, © (3-28) 
and the characteristic impedance is given by 
Z={ 120/Ve(2,+1) ] {Inf (24447 ) /n]} 4 
{25 (7) "inf [2(n41) +4 (a1) 474] /( 2nt4n?iz®) 1 (3-29) 
The results obtained by this method are compared to 
the results obtained in [8], as shown in Figs. 3-2 and 3. 
The method is efficient and the accuracy is reasonable. 
As evident from Fig. 3-2 for €<20 the largest error Sn 


the order of 5%. The accuracy reduces for €=51 to 


approximately 10%, as is evident from Fig. 3-2. For practical 


values of n , the results obtained by Eqs. 3-28 and 29 are 


acceptable. For increased accuracy, more terms in the Taylor 


series expansion of Eq. 3-14 can be included; this procedure 


will be carried out for coupled microstrip lines in the 
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Fig. 3-2 Characteristic impedance of a microstrip 
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following sections. For the sake of comparison, note that 


experimental results obtained by Pucel, et. al. [5a0) indicate 
5 to 10% error for the value of Z over 1 to 6 GHz in the 


widely used Wheeler's formulae fuk 


D. Coupled Microstrip Lines 
Boracouplédimicrostrip lines, Eq. 3-12 and 13° should 


be modified so that the domain of integration includes 


both strips in the coupled lines. Considering Fig. 2-1l0a 


let 
s'=st+(w/2) (3-30a) 
n=w/( 2h), (3-300) 

and 
o=s' /h. | | jee G3-00c) 


Eq. 3-12 can in general be written as 

fo Uy) /iv[1te,/t anh( yh) |} dy=n/2, (3-31a) 
where I is different for the even and odd modes. The 
details are explained for the even modes. The results for 
the odd modes are then summarized. 


For the even modes, 
Iixcos(yx,), aoe) 


where 
Te{ 1/2) S45 Waax( s’ +x) cos( 7x) dx 
+H 1/e) S52 px(s' —x) cos(7x) dx, (3-31c¢) 
and £(x) is the even-mode charge distribution. The charge 
distribution can be expressed as 


p(s' —x) /e=x(s' —x) /M-[2(s' —x) /w]®, (3-32) 


Substituting pf from Eqs. 3=32. into® 3lc,and- introducing «the 
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change of variable 
sind=2(s' —x) /w (3-33) 
to remove the singularities at x=s'tw/2 results in 
Ip=2(w/2) S27 xx( 0. Sw sind) cos[7(s' —0. Sw sin@)]dé, ee 
The slowly varying function xX, can be expanded in powers 
of (s‘'-x) as follows 
(w/2)x(s' —x) =KtK,(s* —x) H(s' —x) #4 -- 
=K,+Kysin@ +K,cos24+K,sin30+--. (3-35) 
Substituting from Eq. 3-35 into 34 and utilizing the 
tabulated integrals of J49] results in 
In=2rrcos(ys' ) Ke m( 0. Sw) 
+27sin(ys' ) & Kerad eva( 0. 57w)- ( 3-36) 
Substituting from Eqs. 3-23 and 36 into 3la, and introducing 
the change of variable 
yvh=a, 
results in 
(e,t1) /2-EF,(-7)™ 
{Re Sof [ee Fe] /ay - 
{+exp[ -j( £s:+¢) a] +exp[ -j( €:-0) a] }Jan(na) da 
+KSslm Jef [ete Fl) 4] /a} 

{ exp[ —j( £40) a] -exp[ +4( s-7) oF] f Jers ( 70) dag Wan) 
where the plus signs hold for the even modes and the 
negative signs for the odd modes. The integral which is 
essential for solving Eq. 3-37 is 

Wn, ¢) =o) (ne) exp[ 4 2n44¢) a] da. 
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The integral can be found in [ 49] and [5°]. For n¥0 and v+#0 


(fE+C) 


ehiswhastdi, Bae .eIF ode ie bie BE a8 ae 


(M-€) 
a 
a 


apis. tueys avben Pe ue vot biked ieute om itt xen 
et do dew kpAeini edt eabioh be ona 1 engin ovitegen. | 
i) pi giao oY taken 

i. abit ys it)~ fem ee ek ae 

OY: bine Catt "0 ee) bie (e+ a et eg al aeetiniel sania 


we ; A 


i ‘A a : yi 
* a : —" 1 ees a fue 
Be Ui a ane ome Ye 


Le 
G 


89 


the result is 
Wn, ¢) =exp[ -jvd(n, ¢)] /{v[ R(n, ¢)]t. (3- 38a) 

R and © are defined as EGINOWS, let 

PCa, ¢) =(2n4i¢) + 2n45g) 247%, (3-38) 
then P(n,@ ) can be written in polar form as 

P(2, ¢) =nR(n, €) exp[jé(n, ¢)]. (3-38c) 
Thus, n-R(n, € ) is the modulus of P, and @ is the phase 
of P. Similar expressions can be obtained for the special 
cases where n=0 or yv=0. Eq. 3-37 thus yields a matrix 
equation in the terms of unknown Ky. From Ky the charge 
distribution and thus the parameters of the coupled lines 


can be evaluated. 


E. Numerical Results 

The results obtained by this method for coupled micro- 
strip lines are presented in Fig. 3-4 and are compared with 
the results obtained in [4]. As in the case of single lines, 
the results are more accurate for small values of w/h (less 
than unity). The charge distribution on the two strips are 
perturbed less for larger values of s/h than in the case of 
Single lines. Therefore, the results are more accurate for 
larger values of s/h. A complete error analysis is 
presented below. 

There are two infinite sums which are truncated in the 
numerical solution, namely, 
Ue the infinite sum of Eq. 3-24 is truncated to N,, terms, 
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Fig. 3-4 Even- and odd-mode characteristic impedance vs. 
w/h for s/h=0.2; coupled microstrip lines with ¢,=9.0; 
Solid curves: [4] ; 
@: present method, even mode, 


@: present method, odd mode. 
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an the infinite sum-of Eq. 3-35 is truncated to Ly, terms; 

Lm 1s the final matrix size. 
The effect of each truncation is considered separately. 
The Effect of Ny 

To assess the effect of Ni on the numerical results, 
La is assigned the fixed value of 5, and Nn is increased 
until the values of the calculated parameters are no longer 
influenced by Nae The results of one typical analysis are 
presented in Fig. 3-5 below. A value of N70 is seen to be 
sufficiently large to produce negligible truncation error; 
in fact, such a value is large enough for all cases of interest. 
thegeitect of Ly 

To assess the effect of Lm on the numerical results, 

Nn is assigned the fixed value of 50, and Ly», is varied. 

The percentage error in each case can be obtained and 

plotted versus Lm. A typical result: is shown in Fig. 3-6. 
Unfortunately, as Lm increases, some of the intermediate 
numbers become excessively large, causing an overflow in 

the computer (specifically R” in Eq. 3-38a). Thus, a 

further increase in the final matrix size is not useful. How- 
ever, for values of less than unity this overflow does 
not occur. 

An example of typical values of the different para- 
meters in this error analysis is as follows: For w/h=1.5 
and s/h=0.3, Nm=30 and L,=9 (the final matrix size) will 
result in 5% error for the values of the characteristic 
impedance. The results will be more accurate for smaller 


values of w/h and larger values of s/h. In general, the 
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Fig. 3-5 


Percent error in final value of Ze vs. Ny; (Ny is 
the number of terms to which the sum of Eq. 3-24 
is truncated). 
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Percent error in’Z) vVs.onl, (Ly, 1s’ the final matrix 
size). 

W/N=le55. S/N=l.0), cn=9.0 and N=40. “(Ny is’ the 
number of terms to which the sum of Eq. 3-24 is - 


truncated). 
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recommended value for L, is 9 to obtain results with 
reasonable accuracy. The method developed in this chapter 
is quite efficient and the computer time and memory space 
required is minimal, since the resulting matrix sizes are 
small (of the order of 8x8). For the sake of comparison, 
the final matrix size obtained by the mode matching method 


is 40x40 as explained in Chap. l, [3]- 


PF. Summary 
In this chapter, an integral equation method was 
developed in order to obtain the potential distribution, ¢, 
at the cross section of microstrip lines in the quasi-TEM 
regime. Green's second identity was considered 
Ss (eV 7y-W7¢) dS= 
Jc [ o( 8y/an) —K ay/an) ] dc; oe 
where. g satisfied Laplace's equation. Unlike Green's 
functions, we also required that 
| vy. | (3-5) 
This reduced the left hand side of Eq. 2-3 to zero. By 
properiy choosing the boundary conditions on Y it was 


possible to reduce theright hand side of Eq. 2-3 to 
"I [p(x) /e] cos (7x) cos( yx) /{71 142,/tanb( yh) ]} dxdy=n/2. (3-10) 
<i 4 


The singularity in the charge distribution was in- 
cluded in Eq. 3-10; | 
p(x) <ex( x) /M{2x/m)* ewe 
P(x) was in turn expanded in terms of a Taylor series with 
unknown coefficients. Thus, the integral on the left hand 


Side of Eq. 3-10 could be evaluated. For microstrip lines, 


the following expressions were obtained for € and 2. 
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ee=f (2,41) /2] Inf (2+447°) /n] / | 
{3 (—7) nl [2( 041) +4(n4t) 247°] /( 24a) Jf, 
Z = 120/r8(2,44) | {inl (244an® ) /n 
{2 (—r) nl [2(n41) +4 (041) 8477] /(2nt4n?in®) 37% (3-29) 


where 


(3-29) 


Sate 1) /( e041), (3-23b) 
The numerical results obtained using the above expres- 


Sions were shown in Figs. 3-3 and 4 and were compared to 
the results of Ref. [8 ]. The accuracy was better than 5% 
for £4 20 and w/h4l. 

For coupled lines, the method led to a matrix equation 
from which Z and €é, could be obtained. The results obtained 
by this method were presented in Fig. 3-4 and were compared 
to the results of Ref. [4]. 

G._ Conclusions 

The method developed in this chapter makes the use of 
Green's functions more flexible. Green's functions are 
obtained by placing a point charge in the geometry of 
interest. The final solution is then obtained by the super- 
position principle. In the method developed in this chapter 
no point charge is needed in obtaining a fundamental solu- 
tion by which the actual solution can be obtained. The 
boundary conditions on the fundamental solution can be 
arbitrarily chosen to simplify the final integral equation. 
The inaccuracy in the obtained results is caused by the 
approximations which are introduced in solving the obtained 
integral equations (specifically, truncating the infinite 


sums of Eqs. 3-24 and 35). Chap. 3 can be viewed as a 
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building block for the wave analysis of microstrips. The 
theory developed for solving Laplace's equation is applied 


to Helmholtz's equation in the following chapters. 
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IV. WAVE ANALYSIS OF MICROSTRIP LINES 


Ae. introduction 

The wave analysis of microstrip lines is presented in 
this chapter. The theory developed in Chap. 3 is applied 
to Helmholtz's equation for the wave analysis of microstrip 
lines. Furthermore, rather than using the two customary 
potential functions ¥° and y™ which were introduced in 
Chap. 1, this chapter uses one function only; it is shown 
that the boundary conditions for EY are separable from the 
other components of E and H. The choice of EY as the 
unknown function greatly reduces the numerical effort of 


analyzing a microstrip line. 


B. Formulation 
The boundaries and the interfaces of microstrips are 
along or perpendicular to the coordinate axes. Therefore, 


the boundary conditions for the electric field, E, and the 


magnetic field, H, are greatly simplified. On the conductors 


Eten=0, (4-1a) 
H..=0, (41d) 
aE,.,/dn=0, (4-1¢) 
oH, .,/dn=0, (4-14) 


where the subscript "tan" denotes the tangential component, 
and "nor" denotes the normal component. 
The continuity of the normal component of the dis- 
placement vector at the DAI implies that 
&2( Eylyas) = Eyl pt’). (42a) 
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On the MIC 
€ (Eyl ya-) + Eyl eat) =p/e, (4~2b) 
In the absence of charge 
V- E=0, (4-3) 
Thus, 
( 8B, /8y| yan-) { OE ,/82+9E,/87] you (4a) 
( 9E,/8y| yan") { OE,/8x49E,/8z] nt (4-trb) 


Subtracting Eq. 4-4b from 4a results in 
( 8E,/8y| yan") -( 9Ey/87|yux:*) = 
(8/8x) [ (Exly-n*) 4 Exlyen”) ] + 8/82) [ (Exlyen*) 4 Eelya) J. (4-5) 
E, and E, are continuous across the DAI, so Eq. 4-5 leads 
to 
( OEY/Oy|y-n-) = SE,/Ay] yon") on the DAI. (4-6) 
According to Eq. 4-lc 
( 3E,/8y|yn-) + 8E,/Oy|,.*) 0 on the MIC. (i= 7 
Eqs. 4-2, 6, and 7 indicate that the boundary conditions 
Lor Ey are separable from the other components of E and H. 


Similarly, it can be proven that 


(Hyly-a-) = Hylya*) on the DAI, (4-8a) 
(Hyly-a7) = Hyjy-n*) =O on the MIC, (4-8) 
( OH,/8y|,-2-) { 8H,/dy|y+) on the DAL (4-8c) 


Eqs. 4-2, 6, 7, and 8 do not prove that pure LSE or 
LSM (Longitudinal Section Magnetic) modes can exist in 
microstrips. Rather, they indicate that microstrips can 
be analyzed by considering EY and Hy separately. For the 
determination of the fields, both EY and Hy should be found. 
However, the quantities of interest, such as &€, of microstrip 


lines or the resonant frequencies of resonators can be found 
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from either Sa or iy and the boundary conditions. The above 
conclusions are based on.the assumptions that the dielectric 
and the conductors are lossless, and that a microstrip line 
is infinitely thin. | 
For the case of microstrip lines, the variation of 
the fields along the z-axis is as exp(-jJBZ ). Thus, Ey 
Should satisfy 

| VeEy+ «k®-6")E,y=0. (4-9) 
In the following sections, the theory developed in Chap. 3 
is applied to solving Eq. 4-9. 


C. An Integral Equation Formulation 


Green's second identity can be rearranged as 
y 


Ssi ¥[ Vigt{ k?-87) 9] -y[ V2yH{ xk?-8?) y] } ds= 


Sel YW 89/8n) —o( ay/Sn) } dc. (4-10) 
where S is the area enclosed by the curve C. Let 
p=Ey, (4-11) 
and let W be an auxiliary function such that 
Vev+{ ck*-B?) y=0. Cia? 


The boundary conditions for w will be described below. 
Considering Eqs. 4-9, 11, and 12, Eq. 4-10 reduces to 
Sel W( 8p/8n) -y( ay/dn) ] dc=0. (4-13) 
The line is assumed to be shielded by a metallic plane at 
y=h+d (Fig. 4-1). In the x-direction, the line is assumed 
to extend to infinity. 
Eq. 4-13 applied to S, and S, in Fig. 4-1 results in 
Scyucgl W 8¢/8n) —p( 8¥/8n) | dc=0, (4-14a) 
Segue,l Wi 9¢/8n) -v( ay/8n) ] dc=0, (4-14) 
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Fig. 4-1 A single microstrip line shielded on top by a 


metallic plane. 
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where it is assumed that @ and a¢@/an tend to zero for 
large x. To simplify the present problem, the boundary 
conditions on Ww are defined as follows: 

8y/an=0 at y=0 and y=hid, (4-15a) 
and 
(8Y/89|y-n-) =er( 89/87 yon*) . / (4-15b) 
Imposing the boundary conditions of Eqs. 4-15 on 14 and 
manipulating yields 
Soa (Hy-n7) 4 Hy-n*) ] (89/8y| yon) dx={ 1/8) Since (x) (99/8 Yl yon*) dx. 
(4-16) 
It is necessary to find ¥ such that Eqs. 4-12 and 15 


are satisfied. Thus let 


rcos(ax)cosh(7' y) /|7' sinh(7' h)]; O<y<h, 
fhe 
os( ax) cosh[y(y-h-d)]/[y7 sinh(yd)];  h<ysht. 

(4-172) 

where 
e,k*—p?+y' ®-o®=0, (4-17b) 

and 

k?-8? +? =0. (4-17¢c) 


Substituting Ww from Eq. 4-17a into 16 results in 
SaarG{ a) cos ( ax) ( 89/8] yan) dx={ 1/8) Sco (x) cos( ax) dx, (4-18a) 


where 
G(a{e,/[ 7’ tanh(7’ h) ]}HY{1/{yt anh( 7d) J} (4-18b) 
Multiplying both sides of Eq. 4-18a by cos(&%xj;) and 


integrating over w from zero to infinity results in 


SoaK( x, x1) (89/3y|y4) dx=0, (4-194) 
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where x;is assumed to lie on the MIC and 


R (x, xi) =2,f_2G( a) [ exp( ja|x+x,|) texp( jalx—x,) ] da. (4-19d) 
Let 
K(x' ) =2nf_2G( a) exp(ja|x' |) da, (4-20a) 
then 
K(x, x,) =K(x+x,) +K(x-x,). (4~20b) 
For a single microstrip line, Eq. 4-19a can be written as 
Sure &( x—w/2) [ K(x+x,) +K( x-x,) ] dx=0, (4-21a) 
where 
a(x—w/2) + 39¢/8y| yon). (4~21b) 
The following changes of variables simplify the above 
equations. 
v=ah, (4-22a) 
¢=x/h, (4-22b) 
n=w/(2h), (4—~22c) 
_K(€) =K(hé), (4-224) 
a(¢) =a(hé), (4-22e) 
G( v) =G( v/h) /h, (laos) 
6=d/h. (4-222) 
Considering Eqs. 4-22, Eq. 4-2la can be written as 
Sova €)[K( €+8,+2n) +K($-€,) ]dé=0; >. ae) 
where (Eq. 4-20a) 
K(£) (2m) f2G{v) exp(jv|é|) dv. (4-23b) 


and (Eqs. 4-17 and 18b) 
G(v) =f 2,/[ ®-C#t anh( W®-c¥) J} 41/1 248% anh(é6W?4B%)].  (4-23c) 
with 
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and 


C®{ e,-2,) (hk) 2. (4~23e) 


D. Evaluation of the Kernel K(é) 


K(&) of Eq. 4-23b can be evaluated by complex inte- 
gration. Assume that y is a complex variable, and integrate 
G(v)exp(jyl@|) along the infinite semicircular path,[ , 
shown in Fig. 4-2; it can be shown that the integrand 
vanishes over x, ; therefore 

K( €) =2nf2G(v) exp(jvlé|) dv= 

-j[(0.5 Res. of poles on the real axis)+{IRes. of poles inside I) ]. 


(4-24) 

The poles of G(v) are located at (Eq. 4-23c) 
po, (4-25a) 
por, (4-25b) 
jog=iMnm) *-C?; n=1,2,3,-°-, (4-25c) 
jo,=jnm/6) *4B", n=O, 1,2,-°°. (4-254) 

Eq. 4-24 leads to 
K( €)e{gfsin( Clé|)/(2C)} 


-{e,5 Jexp( dl él)I/e}{2,lexp( -oal f|V/( 5o)}. (4-26) 
The kernel of Eq. 4-26 is not sufficient for solving 
the integral equation. This is evident by considering the 
limiting case k=0, which implies (Eq. 4-23e) that 
Cheeo=VE rE o( BK) |1-0=0- (4-27) 
For the quasi-TEM mode in a microstrip line, €e remains 


finite for k=O. Thus, from Eq. 4-23d 
Blywo=V6 «1 ( HK) [n=O aes 
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Fig. 4-2 Contour for the complex integration of Eq. 4-23b; 


%, is used to find a complementary kernel. 
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Since (from Eq. 4-25da) 


@g7B-0 as +0, 
the first term in the Berend summation of Eq. 4-26 tends to 
infinity and K(2) becomes independent of C. It is there- 
fore necessary to find a complementary kernel which is 
dependent on C for small C. 

Both sides of Eq. 4-18 can be multiplied by cos( yj) 
and integrated along the small circle around v=C, namely 
3. in Fig. 4-2. vy and ¢ are defined in Eqs. 4-22. 

Thus, Eq. 4-18a becomes 
Sn 1 $,,G( v) cos( vt) cos( vg.) dv] a( £7) dé 
=(1/e) Sue (ht) [ $,,cos(vé) cos(vé,) dv] dé 

=0; regardless of §&,. (4-29) 
The right hand side of Eq. 4-29 is identically zero since 
none of the functions in the integrand have a pole at ys=C. 
Considering Eq. 4-23c G(¥) has a pole at v=C, therefore 
the left hand side of Eq. 4-29 leads to 

So"cos[ C( §+n) ]- a( &) dé=0. (4-30) 

For the quasi-TEM mode, Eqs. 4-23a, 27, 28, and 30, 
prove that €. should remain finite at zero frequency. 
Consider Eqs. 4-23d and 23e. As k tends to zero, both B 
and C tend to zero. The kernels of Eqs. 4-23a and 30 both 
become constant for B=C=0; in this limiting case, Eq. 4-23a 
tends to 

Seral€) dg=0, 
and Eq. 4-30 tends to 
Seval€) dé=0, 
Thus, regardless of the method used to solve Eqs. 4-23a and 
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30, B=C=0 is a solution. This verifies that €e remains 
finite as the frequency tends to zero (from Eqs. 4-23d and 


23e). 


E. Numerical Solution 


In order to solve the above equations numerically, 
f(@) is defined with respect to a(@) of Eq. 4-22e; 
a(£) =exp(-€L)f(€), (4-31) 
where L is an arbitrary number. This definition is based on 
the fact that the fields decay rapidly away from the MIC 
edges. ‘L can be varied arbitrarily to obtain the most 
accurate numerical results; the obtained numerical results 


are nearly independent of L. Then Eq. 4-30 can be 


approximated by 


2, f fefslexp( -€L) cos[ C( ¢+n) ] dé =0, (4-32a) 
where M determines the final matrix size, and 
ft(e). (4~32b) 
Eq. 4-23a can similarly be approximated by 
Bf Sef K(€+ec+2n) +K( €-£) Jexp(—¢L) dg-o, (+32) 
The relationship between W, ,&; ,&;, and: is shown 


in Fig. 4-3; the € axis is divided into segments of width 
W; about the point 4; , ranging from the lower edge éy, 
to the upper edge ¢y. Wj) and &; are the parameters for 
the quadrature integration [ 48]. This choice of parameters 


reduces the integration errors. 


The following definitions will simplify the procedure Tox 


evaluating the integrals of Eqs. 4.32a and 32c. 
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pos | pips 


Fig. 4-3 & -axis divided into various segments for the 


numerical integration of Eqs. 4-23a and 30. 
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Gy, yal 1/(2C) ] Sei sin[ C( §+¢;+27) ] exp(-€L) dé, (4-33 2) 
Ga, yl 1/( 2c) ] ff sin[ Clé-¢,)] exp( -¢L) ae, eon) 
Gan yt 1/pa) Sef! expl ( £+¢:+27) ] exp(-€L) de, (4-33¢) 
Gan u{1/pn) S464 exp[ -palf-€4] exp(-£L) ae, eke 

Gs, =Sefst cos[ C( +m) ] exp(-€L) dé, bee oct 

Goa, y= B/o) Sef) expl -0,( £+£,+2n) ] exp(-€L) dé, Geel 

Gra, y= B/on) Sef! exp[ -o,]¢-¢1] exp(-€L) dé, (4-338) 


The integrals of Eqs. 4-33 can be easily evaluated. Con- 
sidering Eqs. 4-33 we can now rewrite Eq. 4-32a as 

B,f 6s, 0 (4-34a) 
and Eq. 4-32c can be expressed as 

Ets £,OBL Gy gta B (Gye, y+Gan, y) ] “2, ( Gen, ytGrm, y) }=0- (4-34) 

Inwmatrix form Eqs. 4-34 can be represented as 

[G]it}=0, (4-34) 
where the first row of [c] is G,; of Eq. 4-33e. The 
matrix [¢ ] is a function of B and C of Eqs. 4-23d and 23e, 
which are both directly proportional to hk or hf; this 
signifies the fact that, for microstrip lines, the normalized 
frequency fh can be used (i.e., the dispersion relation can 
be plotted vs. fh rather than f). In a similar fashion, 
the strip width can be normalized to w/h as in the quasi- 
TEM regime. | 

The numerical results obtained by this method for 

different representative cases are shown in Figs. 4-4 and 
5 and are compared to the results of [15]. [23] and [52]. 
The results obtained by this method are in good agreement 
with other methods for high frequencies. For low frequencies, 


Ee converges smoothly to the zero frequency limit of the 
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Effective Dielectric Constant, ¢, 


Fig. 4-4 


€&=11.7;  w/h=0.96 


20 40 60 80 
(f-h); [GHz- mm] 


€e of a single microstrip line vs. normalized 
frequency; 

€epis obtained from quasi-TEM methods (8 ] 
and the method of Chap. 2). 

As the frequency tends to zero, the results 
obtained by the present method smoothly 


converge to é€eo. 


109 


100 


a 


te hm ek le ih i gi el 


a 
/ g 


. : _ tan a] : 
} , ; i : 
——E eee Le ee we : | fo 
' " . ssl eb : 


= Coo 


ag ea ne — 


‘pomet aman 5 


vane yy, ‘ a as — P > 5 Pa 
oman Berea +e a ee sien ap PRMD Rr 


ee) ih ae 
jee hh aa 
i of eo 7 : send oat 7 


peeatanane: ev enkt cantata tants # poe , a oa 


[eh litte a2: ta.0u9 aor? “beoihatea Nah” cowie 
VS +Gael Ze boditom ae bom ~ ugha, 7 
ettveer ans jornas od et beanapaignes eg oe o ‘i il 


a 


» Ee 


Effective Dielectric Constant 


Fig. 4-5 


110 


f-h GHz- mm 


€& of a single microstrip line vs. normalized 
frequency; €,=10.1 and w/h=2.5; 

€eo is obtained from quasi-TEM methods (f8 | 
and the method of Chap. 2). 

As the frequency tends to zero, the results 


obtained by the present method smoothly 


converge to Cay. 
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effective dielectric constant, €., , obtained by the quasi- 
TEM methods. The results of [15], [23] and [52] do not 
converge to Cep at zero frequency. 

Considering the numerical efficiency of the present 
method the size of the final eigenvalue matrix used was 
14x14 and fifty terms were needed in Eq. 4-26. The final 


matrix size in the TLM method [22] is 270x270, in the finite 
difference method [19] it is 100x100, and for the finite 


element method 21 it is 126x126. The computer programs are 
presented in App. C. The results of Figs. 4-4 and 5 have been 
obtained by finding the values of & at which the determinant 
of [c] in Eq. 4-34 vanishes. 
F. Summary 
In this chapter, it was shown that, for microstrips, 

the boundary conditions for Ey could be separated from the 
other components of E and H. Ey satisfied Helmholtz‘s 
equation 

ViEy+{ k*—8*) Ey=0. (4-9) 


with the boundary conditions 


SpE Se =n Oe): (W=2e) 
on the DAI, and 
€ ( Ey] ya") {Eyl yen*) =0/E, (4-2b) 
on the MIC. Moreover, 
( 8E,/Ay| yen) = 9Ey/8| y=n*) on the DAI. (4-6) 


and 


( OE #8 isan) =f GE, /OY\ yen*) =0 on the MIC. (4-7) 
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Green's second identity was rearranged as 
Sak ¥f Vigt{ ck?—87) 9] -ol V2y-+{ ck?-8?) y] } ds= 
(Sel W( O¢/8n) -p( ay/an) ] dc. (4-10) 
where 
p=Ey, (4-11) 
The method developed in Chap. 3 was applied to 
this case; it was assumed that 
Viyp-{ ek®-8?) y=0. (4-12) 
Considering Eqs. 4-9, 10 and 12, the left hand side of 
Eq. 4-10 could be set equal to zero. By properly choosing 


the boundary conditions on ¥ , Eq. 4-10 resulted in 


SnasG{ a) cos ( ax) ( 39/89] yan) dx={ 1/2) Sugco(x) cos(ax)dx: (4-18) 


G(a) =e, /[ 7’ tanh(y' h)]+1/[ yt anh( yd) J, (4-18) 
where 

Sree ae: (4-1Ta) 

KBP +7270. (4-17 b) 


Eq. 4-18a multiplied by cos(@-x;) and integrated over « 


in the complex a-domain resulted in 


Sore( €) cos[ C( +n) ] dé, (4-30) 
and 
Sonal €) (KC E+€, +27) +K( €-€,) ] dE =O, (4-23a) 
where a(@) was defined by 9 @(x,y)/9 y on the DAI, where 
&=x/h | (4~22b) 


and the kernel was given by 
K( £) =e,sin(C|é|) /(2C) 
—e,5, exp( palé|) /ea-EexP(—oalél)/(don)- (4-26) 
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B=\é,-1(hk), 
C=\é,-e,(hk). (4-23e) 


pam n7)®-c4, n=1.2,9.°"". (4-25c) 
o,—Knr/8) 94B2, NO. 12: (4-25d) 


Eqs. 4-23a and 30 yielded an eigenvalue matrix equation 


(4-23d) 


when discretized. The numerical results obtained by this 
method were presented in Figs. 4-4 and 5 and were compared 


to the results of foie 123] and [52]- 


G. Conclusions 

It has been shown that the wave analysis of microstrips 
can be performed using only E.. iin contrast to the customary 
two potential functions ple) and p(t) the use of only 
one function greatly reduces the effort for analyzing micro- 
strips. In order to solve Helmholtz's equation for Ey 
the theory developed in Chap. 3 was used. The results 
obtained for microstrip lines show good agreement with 
other methods (Figs. 4-4 and 5). As the frequency tends to 
zero, the results converge smoothly to the values obtained 
by quasi-TEM methods. The method is economical (final matrix 
size of 14x14). However, it is not readily apparent how 
this method can be generalized to microstrip resonators. 
Therefore, in Chap. 5 the method of this chapter is modified 
so that it can be generalized more easily to microstrip 


resonators. 
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V. AN ALTERNATIVE INTEGRAL EQUATION METHOD FOR 
THE WAVE ANALYSIS OF MICROSTRIPS 


A. Introduction 

It seems difficult to generalize the integral Slersina 
method of Chap. 4 to microstrip resonators of arbitrary 
shapes. Considering Eq. 4-19a, the domain of integration 
of this integral equation is the DAI. The resulting matrix 
equations would be smaller if an integral equation could be 
found whose domain of integration is the MIC alone. The 
generalization of such an integral equation to microstrip 
resonators would then be simpler. The development and solu- 
tion of an integral equation for microstrip lines with only 
the MIC as the domain of integration is discussed in this 
chapter. This solution is then generalized to a rectangular 


resonator, resonating in the even mode. 


be hoOrmiuLation 


Eq. 4-18b can be inverted to result in 


1/G{ a) =7- sinh(7d) 7’ - sinh(7' h) / 
[e,7- sinh( 7d) cosh(7' h) +7’: sinh(7’ h)cosh(yd)], (5-1a) 
where 
y=" +p7—-k*, (5-1b) 
of Boat +p?—e pe’, (5-1c) 
Multiplying Eq. 4-18a by Eq. 5-la results in 
Sanyoos ( ax) (89/8y| yu) dx={ 1/2) Syacl cos (ox) /G(a)]p(x)dx. (5-2) 
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where es Eq. 5-2 cannot be integrated over ¢ from 
zero to infinity, since 
1/G(a)~a for large a 

However, both sides of Eq. 5-2 can be multiplied by 

cos(ax,)- sin*( ad /2) /a?, 
and integrated over a from 0 tom, 
- Thus, the left hand side of Eq. 5-2 results in 


Sony £08) 2/aY |.) ax (5-3a) 


where 
f (x) ={)"cos( ax) cos( ax,) sin®( ad’ /2) /a®da (5-30) 
f(x) of Eq. 5-3b is nonzero only on the regions 
x,;-d, <x<x,+d, , and 
(xytdy ) <x<{x,-d ), 
shown in Fig. 5-1. Thus, f(x)=0 on the DAI if 
|x,+d, |<w/2 V i. (5-4a) 
With the condition of inequality 5-4a, the integrand of 
Eq. 5-3a is identically zero. Thus, the right hand side of 
Eq. 5-2 results in 
SucK( x, xy) p(x) dx=0; (5-4b) 
where 
K(x, x) =f9"[ sin?( ad’ /2) /a®] cos( ax, cos( ax) /G( a) da; (5-4e) 
The integral which is needed for evaluating K(x,x, ) 
of Eq. 5-3b is 
K,( ¢) =fo"exp( jv) /[ v*G( v) ] dv. (5-5a) 
where 
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Fig. 5-1 Microstrip line shielded on top by a metallic 


plane; f(x) of Eq. 5-4 is nonzero on -(x,+d,) to 


-(x;-d,), and (x;-d';) to (x, +d). 
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&=x/h; 
In other words, 


6 
K(x) =2 Bq Ki (Xai /h) , 


(5-5c) 


(5-5) 


where the coefficients 4q and x,, can be found easily from 


Eqe 5-3b. Moreover 


G( v) =€( v/h) /h. 


G( ¥) can be rewritten as 
G(v) =Pa(p) /Pi(e). 
where, from Eq. 5-la, 
P,(p) =-psin(p) vA?-p*sinh( 6A?) ; 
Po( p) =e, VA®-p*s inh( 6 VA?—-p*) cos(p) -psin(p) cosh( 6-vA?-p®) 
with 
6=d/h, 
p*={7' h)* 
A?.={ e,—1) (hk) ®. 
pact". 
ch f,—4) (hk) 5 


To evaluate K, of Eq. 5-5a, it is assumed that v isa 


(5-5e) 


(5-6a) 


(5-6b) 
(5-6c) 


(5-7a) 
(5-7b) 
(5-7c) 
(5-74) 
(5-7e) 


complex variable and the integrand of Eq. 5-5a is integrated 


over the contour shown in Fig. 5-2. The integral vanishes 


over the semicircle, Cy. » when the radius of the semicircle 


becomes infinitely large. Then 
K,(¢ ) =nj2[ (Res. on the real axis) 


+{Res. on the imaginary axis) ]. 


(5-8) 


The real poles in Eq. 5-8 are treated differently from 


those in Eq. 4-24 of the previous chapter. The reason for 
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Pig. 5-2 Contour for the 1ntegral.ot Eq. 5-5a- 
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including the real poles inside the contour of integration in 
Fig. 5-2 is that f(x) of Eq. 5-4 is defined on the basis of 
step functions; in defining the step functions by Fourier 
transform methods the contour of integration approaches 
the real axis from below and includes any pole on the real 
AX1S 

It is then necessary to find the poles of the integrand 
in Eq. 5-5a; they can be obtained from the roots of P, (/) 
of Eq.,5-6c. It.can be, shown that, all.the.roots of Eq. 5-6c 
are real. Let the m-th positive root of Py (fP) be denoted 
by ae The poles, Yn? can be found from yi via Eq. 5-7d. 


The poles of the integrand of Eq. 5-5a are located at 
v,=0, double pole; and (5-9a) 


Veg 2p 2, ( 5-9b) 


or 
$y ,=tVer—p ft for py<C; m=1,°-°,M 
Vmn= 
Vm=jp2—C® for p_>C; m=MH,-->. (5-9¢) 
Eq. 5-8 leads to 
K,(¢ )=Q.¢ + 
ey 2Q,sin( vf ) +3 Qnexp(—Yme )- (5-10a) 
where 
Qo={ 17/2) P,(C) /P2(C) - (5-10b) 
Qn=7P (Pm) /L VanPs( Pm) J - (5-10c) 
Qn=-TP (Pn) /L VPs Pm) J - (5-104) 
and 


Ps(p)={1/p) Pel (p) suf Be) 
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C. Single Microstrip Lines 

Due to the even symmetry of the quasi-TEM mode in 

microstrip lines, Eq. 5-3a can be rewritten as 

S0"/*p( x) K(x, x,) dx=0; 

O<x,+d,' <w/2. (5-11) 

where K(x,x;) is a linear combination of forms like K,(é). 
As evident from Eq. 5-10a 

K,(¢€)~in(€) for small ¢; 
K(x,x;) behaves logarithmically as x approaches 

(xd! ), x, or (x;+d) ). 

Moreover, as discussed in the previous chapters, p(x) is 
singular for x=w/2. 

By devising a numerical method which includes both the 
Singularities in the kernel and the charge distribution of 
Eq. 5-11, accurate results could be obtained. However, 
generalizing such a method to microstrip resonators of 
arbitrary shapes seems very difficult. Two approximate 
methods for solving Eq. 5-ll can be devised. The singularity 
an the’ charge distribution ean ve sinrcluded in the kernel’ ov 
Eq. 5-11 and the logarithmic singularities of K(x,x;) can 
be "avoided" by dividing the domain of integration into 
intervals which have | 

(xpd), x, or (x;+d,' ) 
as their end points. 

Alternatively, the singularity in the charge distribu- 
tion can be implicitly included in the integration procedure 


by choosing finer intervals near the singular point, x=w/2. 
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The two methods yielded very similar numerical results and 
the computational efficiencies were also comparable. However, 
for reasons of conciseness, only the former method will be 
explained. 
As in the previous chapters, let 


p(x) =x(x) /M{ 2x/w)*, (5-12) 


Substituting Eq. 5-12, together with the change of variable 


sin@ =2x/w, (5-13) 
into Eq. 5-11, results in 
So"/®x( 0. 5w sin@) K( 0. 5w sind, 0. 5w sind,) dé=0. (5-14) 
The change of variable of Eq. 5-13 removes the singularity 
ate xew/ 2. 
Eq. 5-14 can be approximated by 
¥ xrSafK( 0. Sw sing, 0. 5w sind, d6~0, (5-15a) 
where 
Xr=X(0. 5w sind,); (5-15d) 
sind ,=2(x,—d,' ) /w; (5-15c¢) 
sin@,=2x,/W; (5-154) 
sin@,,=2(x,+d,' ) /w. (5-15e) 
Eqs. 5-15c, d, and e are graphically represented in 
Fig. 5-3. 
Eq. 5-15a can be written in matrix form as 


CH] ix} =0, (5-16) 
where [H] is a square matrix, with elements 
H, =Jof7K( 0. Sw sing, x, dé, ool) 
and {X} is a column matrix with X, as its elements. 


Substituting from Eq. 5-5d into 17 results in 
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Fig. 5-3 MIC divided into several intervals for the 


numerical integration of Eq. 5-15. 
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HB, agfog™Ke( fq) 0 (5-18) 
From Eqs. 5-30, Jd, and 13 x,;; can be found as 
X1,/(w/2) =|sin6-sin6,); (5-19a) 
¥q/(w/2) 4sin6-sing|; (5-19b) 
Xg/(w/2) 4sind-sindy: (5-19c) 
Xu/(w/2) =sind-sind,|: (5-19d) 
X51/( w/2) +sin6-sind|: (5-19e) 
Xe/(w/2) 4siné-sing a); (5-19f) 


Finally K,(2) must be evaluated from Eq. 5-l0a. To 
simplify the evaluation procedure, the following definitions 


are introduced. 


nlst+s’| for m=0, 


f,m(s,8’)=4 2sin(nvygists'’|) for 1<m<M, (5-20a) 
exp(-7v,,/9+8'|) for m>M 
fo m( 8.8") fy (8, 8") a(S, 8"): (5-20b) 
fs m(Siné) +f, (sind, sin®y) 
—2f s (sind, sind, +f m( sind, sind); (5-20c) 


Substituting from Eqs. 5-10a and 20 into 18 results in 

Hip Bs, Qan’ [Soff a, ml SinGO +S of 5, mu(sind) dO] . (522%) 
The integration domain of Eq. 5-21 is divided into two 
regions to avoid the logarithmic singularity in K,(é), as 
discuteet above. The integrals of Eq. 5-21 can be carried 
out numerically using the Gauss quadrature formulae. The 
programming details are presented in App. D. The dispersion 
relation is obtained by the determinant search of [H] in 


the eigenvalue Eq. 5-16, 
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Having formulated the above modified method for analyzing 
microstrip lines, we now generalize it to microstrip 
-resonators. The numerical results for both cases are 


subsequently presented. 


D. Microstrip Resonators 


In this section, an example is given of how the general 
problem of microstrip resonators can be analyzed. The case 
of a rectangular resonator (Fig. 5-4) with the even mode of 
resonance is studied since, for it, the results of the 
previous section can be more easily generalized. Eq. 5-3a 
can be generalized to 

SgsK( x, 2: &y 2) (x, z) dxdz=0 (5-22) 
where ((x,zZ) is now the two-dimensional charge distribution 
on the MIC and 

K (x ,Z3%p2i) =SoSo ( P1/P2) (sinad,’ /a) *( singd," /f) *. 
cos( ax) cos(ax,) cos(8z) cos( Bz,) dadg. (5-23) 
P,(P) and P(P) are defined in Eqs. 5-6b and 6c. From 
Eqs. 5-le and 5=-7b 
p®=D*~ ah) °{ gh) *, (5-24a) 
where 
D*=e,( hk) ®, (5-24d) 
K of Eq. 5-23 is a linear combination of terms like 

Kal € ,¢)=S 2/4 Pi/P2) exp[j(vé. +v' ¢ ))/(vv' )*®dvar' , (5-25a) 

where 
€=x/h; (5-25b) 
¢=2z/h; (5-25c) 
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Fig. 5-4 Rectangular microstrip resonator with even mode 


of resonance; + indicates the charge polarity. 
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v=ah, 
(5-25d) 
v' =6h, (5-25e) 
To integrate Ko» the following transformations are 
introduced 
v=ocos6, ( 5~26a) 
v' =osiné, (5-26b) 


Substituting from Eqs. 5-26 into 25a results in 


Ke SXPi/Pz\eyo) do. (Se27e9 
where 
1(6)=2fo*"exp(joE) /[ (1-cos6' )(1+c0s6')]- dé’, (5-27b) 
with 
Q° =20, (5-27c) 
and | 
E=€ cos(0' /2)+{ sin(6' /2). (5-274) 
Eq. 5-27b can be evaluated by complex integration methods. Let 
Z=exp(jé' ), (5=28) 
then 
aZ/(jZ) =d0', (5-29) 
and 
I={ -8/}) $eZ- exp( jo) /[(Z-1)(Z+1)]*aZ, (5-30) 


where C is the contour shown in Fig. 5-5. There are two 
double poles at 

Z=1 and Z=-1. 
Whether these poles should be fully enclosed within C of 
Fig. 5-5 or not is a valid question. However, in either 
case, the solutions to Eq. 5-3a will be the same, since the 
right hand side of Eq. 5-3a is zero and multiplying both 


sides by a constant will not affect the solutions. 
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Fig. 5-5 Contour for the integration of Eq. 5-30. 


128 
To evaluate I of Eq. 5-30, the residues of the integrand at 
the right hand side of Eq. 5-30 arenes found: Let 
Ry 4/42) [Z- oxp(JoE) /(Z+1) "Vay 
={[ exp(joE) +Jo( dE/do' )( de" /aZ) Z- exp(joE) ](Z+1) * 
—2(Z+1) Z- exp(joE) } /(Z+1) ‘|z.1. (5-31) 
From Eq. 5-28 
6" |g_4=0. (5-32) 
substituting from Eq. 5-32 into 27 results in 


E|za.=0. 5€, (5-33a) 
and 
dE/d6' |z4=0. 5¢, (5-33b) 
substituting from Eqs. 5-33 into 31 results in 
Ry={ 1/8) of exp(joé.). so) 
Similarly, 


Re={ d/dZ) [ Z- exp(joE) /(Z-1)*] z-4={ -1/8) o€ exp(jog). (5-35) 


The addition of the two residues yields 


I= -8/§) (27) (Ri +Re) 
=16no[¢ exp(jof )+€ exp(jog )]. (5-36) 
Substituting this result into Eq. 5-27a, we obtain 
Ke=<M( £) -€M{ ¢) (5-378) 
where 
M( £) =/(P1/Pz) [exp( jot) /o*] do, (5-370) 


The evaluation of M(&é) is similar to the evaluation of 
Kj (¢) in Eq. 5-5a, above. The results are similar to 
Eq. 5-10. In this way, Ko Or WEGs. 5-57a-1s found, 
The kernel of Eq. 5-23 is a linear combination of terms 
like Ky of Eq. 5-37a. ‘Having Tound* the kernel, Eq. ©5=23 


can be discretized. The resulting eigenvalue matrix equation 
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yields the resonant frequency of the resonator. The procedure 
for finding and solving this matrix equation closely parallels 
that in the previous section. The numerical results obtained 


by this method are presented below. 


E. Numerical Results 

The numerical results obtained by this method for 
Single microstrip lines are shown in Fig. 5-6 and are 
compared to results obtained in [24] . The results show 
good agreement with [24] for higher frequencies. This 
method is less accurate for lower frequencies. At f.h=20 
GHz-mm, for example, the error is 10.7%. 

A numerical error analysis is performed as follows. 

The effect of truncating the infinite sum of Eq. 5-10a to 
N, terms is demonstrated in Fig. 5-7. The final matrix 
size, Mz, is assigned the constant value of 4. The curve 
Derk the change in the calculated value of ¢€ vs. Nx; 
the kernel of Eq. 5-11 converges very slowly. Fig. 5-8 
shows the effect of increasing the final matrix size, M,. 
N, is assigned the fixed value of 500. It is seen that an 
increase in the matrix size from 4x4 to 8x8 has little 
effect upon the obtained results. 

Unfortunately, the rate of convergence varies for 
different frequencies. For frequencies below the inflection 
point in the dispersion curve, even values of Nx,=1000 provide 
poor results. Furthermore, the cost of running the computer 
programs becomes prohibitive as Nx becomes very large, even 


though the matrix size is as small as 4x4. 
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Fig. 5-6 Effective dielectric constant of a microstrip line 


vs. normalized frequency; ¢,=11.7 and w/h=0.96 
Solid curve: method of Chap. 4, 
Dashed curve: Ref. [24], 


Plus signs: present method. 
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Calculated Effective Dielectric Constant, ¢, 


tO. O 
9. 8 
9. 6 
9. 4 
Fig. 5-7 
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400 © ©200.—~2~=—*«é‘é SCO 400 


Calculated effective dielectric constant 
vs. Ny (Ny is the number of terms to which 
Eq. 5-10a is truncated). 

fh=40.0 GHz mm, w/h=0.96, &=11.7 and 
Mp=4 (Mp, is the final matrix size). 


The actual value of & is 10.5. 
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Final Matrix Size, Mp 


Calculated effective dielectric constant vs. 
final matrix size, Mg; fh=40.0 GHz-mm, w/h=0.96, 
& =11.7 and N,=500. (NL is the number of terms to 


which Eq. 5-10a is truncated), 
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Thus, the approximate range of usefulness of this 
method is the frequency range above the inflection point 
in the dispersion curve (Pig. 5-6). In order to obtain 
accurate results with this method, the dimensions of a 
rectangular resonator should be chosen such that the 
resonant frequency remains within this range. In other 
words, the product of the resonator's expected resonant 
frequency and substrate height should be above the 
inflection point in the corresponding dispersion curve 
(i Jes, the disperaion curve for the line which results 
from extending the resonator length to infinity). 

Thus, a resonator with the following specifications 


is analyzed (Fig. 5-4): 


é,=11. 7, 
h=3. 17 mm, 
w/h=0. 96, and 
L/h=3. 263. 


According to the transmission line model described in Chap. 
1, where dispersion in the line is included and fringing at 
the ends neglected, the resonant frequency should be 9.2 GHz. 
With the method developed in Chap. 5, the calculated 

- resonant length for this resonant frequency is L/h=3.248. 
This indicates an error of 0.15% from the specification. 
However, the transmission line model by itself may be in 


error by a few percent [32] and [35] . 
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EF. Summary 
In this chapter, the integral equation method of Chap. 
4 was modified such that the domain of integration of 
the integral equation was the MIC. In this way, the method 
could be more easily applied to the rectangular resonator. 
Eq. 4-25b was rearranged as 
Soarcos ( ax) (89/3! yn) dx={ 1/2) Sacl cos (ax) /G(a)]p(x) dx. (5-2) 


where ¢=E,and 


1/G( a) =y: sinh( yd) 7’ - sinh(' h) / 


[ ¢,y: sinh( yd) cosh(7' h) +7’ - sinh(y'h)cosh(ya)], (972?) 


with 


y=a* +p? 18, (5-1b) 
7’ *=08+p*—e ke", (5-1c) 


Both sides of Eq. 5-2 were multiplied by 
Cc - sin? g 4 
os ( ax,) sin*( ad‘ /2) /a?, 


and integrated over a from zero to infinity to result in 


SacK( x, x,) p(x) dx=0; ( 5-3a) 
K(x,x; ) was a linear combination of forms like 
Ky( ¢) =fe"exp( ive) /L v76(v) ] av. ee 
where 
opie (5-5) 
(=x/h; (5-5e) 
Gv) =G{v/h) /b. (5-5a) 


K, (&) could be found by complex integration techniques. 


Eq. 5-3a was further developed for the case of the 
rectangular resonator in the even mode of resonance. 


The numerical results for microstrip lines obtained by 
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this method were presented in Fig. 5-6 and showed good agree- 
ment with the results of [24] at higher frequencies. The 
method failed to be accurate at lower frequencies. The 
range of usefulness of this method was for frequencies 
higher than the inflection point in the dispersion curve. 
Using this range of usefulness, the method was then applied 
to the rectangular resonator and a numerical example was 
provided, which showed good agreement with transmission line 


model methods. 


G. Conclusions 

The method developed in this chapter reduces the problem 
of analyzing microstrips to an integral equation with only 
the MIC as the domain of integration. The method for solving 
the final integral equation has been presented. This method 
gave results which are in good agreement with those of other 
methods for high frequencies, but is less accurate for low 
frequencies (e.g., the error at f-h=20 GHz. mm is 10.7% 
for the dispersion curve shown in Fig. 5-6). 

K,(¢) of Eq. 5-10a converges very slowly and becomes 
logarithmically singular for small arguments. The numerical 
analysis indicates that a better accuracy can be achieved 
for lower frequencies if more terms are included in Eq. 5-10a. 
However, this causes a large increase in computational costs 
even when the final matrix size is as small as 4x4. For 
frequencies above the inflection point in the dispersion 
curve, the method is quite useful and can be generalized to 
microstrip resonators, as indicated by the example given in 


this chapter. 
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VI. CONCLUSIONS 


In the previous chapters, various integral equation 
methods were applied to the analysis of microstrips. The. 
BEM was applied to single and coupled microstrip lines in 
the quasi-TEM regime. The analysis was capable of evaluating 
the characteristic impedance, Z, the effective dielectric 
constant, €g , as well as the charge density, P(x), and the 
potential distribution, g. The inclusion of the singularity 
of the charge distribution in the formulation resulted in 
increased efficiency, increased accuracy, and reduced final 
matrix sizes. 

However, complications were evident for generalizing 
the BEM to the wave analysis of microstrips, since it was 
difficult to obtain the "fundamental solution" (see p. 45). 
It was easier to obtain the results by modifying the classical 
integral equation solution of Chap. 2. Thus, in Chap. 3, the 
quasi-TEM microstrip line was analyzed by a novel method which 
can be considered as a generalization of Green's functions. 

Green's functions are customarily defined so that they 
satisfy a nonhomogeneous differential equation. For example, 
the forcing function can be chosen to be a point charge. In 
Chap. 3, it was shown that instead of Green's functions, a 
fundamental solution could be chosen which satisfies a homo- 
geneous differential equation. The boundary conditions on 
this fundamental solution could be arbitrarily chosen in 


such a manner that the final integral equation was simplified. 
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This method resulted in reasonably accurate expressions for 

the characteristic impedance and & of ahensen aie lines. 

This method was then generalized in Chap. 4 to the wave 
analysis of microstrip lines. The problem of analyzing 
microstrip lines reduced to solving an integral equation 
with the DAI as the domain of integration. Moreover, 
unlike the customary employment of two potential functions 
for the wave analysis of microstrip lines, it was shown that 
only one properly chosen function, namely Eye was sufficient 
for obtaining the parameters of the microstrips. This 
approach reduced the effort for analyzing microstrips. 

The obtained results exhibited good agreement with 
other methods and converged smoothly to the quasi-TEM 
results. Furthermore, there were two conclusions of interest 
regarding the wave analysis of microstrips: 

1. §as in the quasi-TEM regime, the dimensions of micro- 
strips could be normalized to the substrate height, h, 
and 

2 the frequency, f, could be normalized to f-h. 

In order to generalize the formulation to resonators, 
an alternative integral equation method was developed in Chap. 
5 with the MIC as the domain of integration. This method 
was more elaborate, but could be more easily generalized to 
microstrip resonators. The obtained results were accurate 
for higher frequencies but the accuracy was reduced for lower 
frequencies (i.e., below the inflection point in the dispersion 
curve). 


The rectangular microstrip resonator with an even mode 
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of resonance was then analyzed using this method. This 
analysis indicated how the method could be generalized to 
microstrip resonators in order to provide a basis for the 
computer-aided design of microstrips. For the further 
development of the above methods, the following recommenda- 
tions can be employed: 

The general problem of analyzing microstrip resonators 
can be approached in a way similar to the analysis of the 
rectangular resonator. The MIC area of an arbitrarily 
shaped resonator can be subdivided into a mesh of rectangles; 
in this way, the obtained integral equations can be discretized. 
Due to the singularity in the charge distribution, smaller 
rectangles can be chosen near the MIC edges. Furthermore, 
the numerical errors caused by the slow convergence of the 
obtained functions could perhaps be reduced by expressing 
these functions in terms of more rapidly convergent functions 
(e.g., a logarithmic function added to a Taylor series). 

The methods developed in this thesis can also be 
used for analyzing coplanar microwave integrated 
circuits (i.e. circuits in which the dielectric substrate 
is metalized on one side only). Another possible field 
of application can be the analysis of the higher order 
modes in  wmicrostrips. 

This thesis can therefore be viewed as a stepping stone 
toward the eventual analysis of arbitrarily shaped resonators. 
In addition, the methods proposed for the analysis of micro- 
strip lines and coupled lines can be applied directly to 


engineering designs. 
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APPENDICES 


Appendix A 


Computer Programs for Single MIC Lines 


The programs are written in APL. The workspace is 


called SINGLE. The symbols used in the programs and their 


relation to the symbols in Chap. 2 is as follows: 


CH 
CO: 
Ci 
ER? 


{X} in Eqs 2-18 and 37 

Co capacitance of the microstrip line for f€y=l 
C capacitance of the microstrip line for €71 
er 
4 Eq 2-2la 

[c,] Eqs 2-22b and 37 

[G5] Eqs 2-22b and 37 

[G3] Eq 2-37 

[Gy] Eq 2-37 

[G.1 Eq 2-37 

¥; Eq 2-25a 

¥; Eq 2-36d 

Abscissa for Gaussian integration 
@ Eq 2-29c 

6, Eq 2-29a 

6 Eq 2-29b 

W; Eq 2-24 

W, Eq 2-36c 


Weight factor for Gaussian integration 
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WOH: w/h 

Lee G6 Witne O<G < | 
Ei wc. Wiehe scl <e 

ZEFs C* with | 0<¢ 


Program POINT 


This program assigns the values to the parameters 


necessary for numerical integration; the data are obtained 


from Ref. \48]. 


v FOINT 
ba NJ¢ 0.07652652109999999 0.2277858511 0.3737060887 D0 
we 0.51086770019 0.6360536807 
C3 NJ¢NJ+ 0-7463319064 0.8391169718 0.91223544282 5D 
C4] OV GOS27AI 2/2 
Ca NJI¢NI20.9931285991 
C6] WJe 0.1527533871 0.1491729864 0.1420961093 0.1316884384 
EZ WJeWt> O-11TS1945319 O8101SS01198 0.0832767415 
Ce8I WJetWJs 0.06267204830000001 0.0406014298 0.0176140071 
Eo) THLE00.5x0+0.5X CL VNIJDF14N5I 
Ci0J TH2¢C1i¥TH1)+00.5 
C113 THIeEoO.SxXNJ 
Ef21 ZeerotHul 
ERS] Gheéte4 
C142 NFe 0.26354603197 1.4134030591 3.596425771 7.0858100058 +0 
Eis) - «6408008442 
C16] WFe 0.6790940422 1.6384878736 2.7694432424 4.3156569009 3D 
bl71) Ze 2191863544 
Mio Ar GN 
bi 
Program INT 
This program evaluates the elements of [G51 (Eq. 2-32). 
vy RREINT 
£13 LtHe THe bJIS THe 
Geel THEO. S*X:GhHSEI ISTHE) +NMXOTH 
i352) AtTH-THICII 
C4] REO. SXTHtTHICII 
Gol CtB-00.5 
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C6] AR¢€@((10A)=A)X(C10OB)+B)X(10C)+C 
C7 J DL¢DTHx+/WMxAR 
C8) ALECTHICJI-THICII 
C9I a2elnelJI- THe la 
Brod <LeTHLCJI+THIET I 
ClidJ B2¢TH2CJI+THICII 
Ci2] Cit¢ki-o1 
Ci33] C2¢B2-o01 
C14} D2e2xe((CIAZ)*«AZ)XCCEBR2)4AB2)xX0¢102) 402) 
CiSJ DB2en2-2xe(ClAL)*AL)XCCIBID ABI) XC ICI) «Ci 
C16) RReD1+D2-0THXCF 
v 
Program INT2 
This program evaluates the elements of 165) (Bdie2=37). 
v RS¢INT2 
Cid PEG PH2CJI-THIEU 
E22] THEO. SX (CTH2DLIJI+THIC II) +NMXDOTH 
Es) STHtE1OTH 
C4) R1Ii¢STH=STH+ZFCII 
C5] R2¢€¢(CTH-00.5) x20TH)+ZFCII-STH 
£6] RS¢2xXTTHX+/WMXRIXK2 
Le The El aS) al Ga FS A 
eet R5¢(ZS-CLOTHICIJI) 42) a4 THILII-900.5 
co) R6¢€(ZS-CLOTH20C II) 42) 4 TH20I1-900-5 
C102 R7€2xXeRS6=KS 
C11] RS¢R3+R7 


Vv 


Program FIRST 


C1) 
C2) 
C3] 
C44 
C53 
C6) 


This program forms the matrix {¢,1 (Eqs. 2-24, 22b). 


vFIRSTCOIv 


7 FIKST ET 
CEleCZzE TD #2 
JMe PW 


ARI ¢CF1+¢(ZEe.+ZE)k2 
AR2CCEI + (ZE ca ZED KS 
ARS+¢@ARIXAR2 
G1¢-00.S5XARSX CIM» IM) PW 


Vv 


~ 


COMRSewema) cae (aN amesoge 
ba O42) bcd Rae TMP eL AAC LAs Pee 


eat ‘Sees io 
OK 


/ tf 


(.O] to siremeis page ast au Lui wah af aint 


4 fits 


CL DEC ISNT OHTA 
Tuer! ti, sesivihlbdetinnaidh 


| . cry pers ase 
AT2- cow geen 
ohm LAN + KTS 


Senta 
24,00~ be {LTA CARLOS BEER 
¥ See erenereen {TOLD -BSPAR 


s* u ** &. 
* ie) 


»(dSs. .¢§-s sept) f(y veniet ote. oiradt wae (out 


Program SECOND 


This program forms \G,\ using INT2 (Eqs. 2=92,..220). 


v SCNI 
Gt) G2¢(JMs IM) P0 
C2] I¢i 
C32 L2iJel 
C49 L3zG2CIsJICINT 
CSI Je J+1 
C6) 3( J£IM)/LS 
C7I ITeI+1 
C8) 9¢(TSJM)/L2 
Vv 
Program THRD 
This program forms 1G.) using INT2 (Eqs. 2-37). 
vy THRE 
Et] JFep ZF 
e3 G3¢(CJUF »>JM) PO 
Tel 
C4) L4iJei 
ESa Lae GS3ETs JJeINT2 
C6] pel Kg Baa 
C73 2(JZ£IMI/LS 
ce FeLi 
Loe 28 4(TH£JF)/L4 
Vv 
Program FRTH 
This program forms |G, (Eqs. 2-24, 37). 
vOPNRHO ET 
Cig Ch1l¢¢2=-ET) «2 
bed RIGCChIACZE eo. tZb)x2 
Ba ReCGR T+ CZ he. 2 bk ke 
C4) R3¢@R1XR2 
Bsa G4¢-00.5XR3X (UF + JM) PWS 


Vv 
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Program. FET 
This program forms 14.) (Bas. 2=36, 37). 


v FFTH 
bid ZEFCZE > ZF 
C2) JN€¢ Pp ZEF 
C33 R1l¢€+CF1+(¢(-ZEF)°¢.+ZF) «2 
C4] R2¢+CF1+(ZEFo.,+ZF) a2 
CSI R3¢R1+K2 
C6] GS C2xCER-LIFETIXRSX CIN» JF) OWE 


v 


Progran FNL 


This program performs the final assembly for the case 


€.=1 
and solves the resulting equations to obtain {X} and C, 


(ani 2-3h),, 


vy FNL 
Cig GG¢G1i+G2 
C27 FFe JMP04 
£37] CHEF F EGG 
C4) COe+/CHxXWJ 


Vv 


Program FNL2 


This program performs the final assembly for the case 


é,¢1 
and solves the resulting equations (Eqs. 2-37,38). 


v FNL2 ER 
Eid CF SeO2xER+1 
C2] ito CE St CER 1x SOR TX CL-ZER)S2)4+ SOE TX CIFZEF )=2 
C3 Ge((G1i+G2),C011 G3+G64)»G5 
C4] Te JM+1 
CSI Lé6:GCIsII¢GlIslI+CF3 
C6] QPELICFOCII+CES 
C7 Tei+l 
C8] 9CTZIN)/L6 
C9 CHF ¢F RG 
C103 Cie+/( IM*CHF) XW 


v 
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Program TTL 
This program produces the titles needed for the output 


in the main program IMPD. 


beam 
£id Ti¢ ‘CHARACTERISTIC IMFEDANCE' 
C2) Vee CPPECTIVETINIELECTRIC sCONSTANT * 
C3) T3¢ ‘CHARGE DISTRIBUTIONsCH VS. ZE' 3 
C4] T4¢*FPOTENTIAL DISTRIBUTION ON THE IAT! 
CS4 TS¢ FIRST COLUMN: NORMALIZED DISTANCE' 
C46] T6é¢'SECONDI! COLUMN: FOTENTIAL' 
CAI T7e*! 
C8) T8¢'*SECOND COLUMN: CH' 


v 


Program IMPD 


This program is the main program; it executes the above 
programs and prints out ¢€,, 2, charge distribution (X of 
Eqs. 2-18 and 2le), and the potential distribution on the 
DAI. The input to IMPD is w/h and ¢€, . 

vy WOH IMFD ER 


i a ET¢&WOH=2 
C2) FIRST ET 


C3 FNL 

C4) ERT «Ed, 
C34 FFTH 
C6] FNi2 ER 
C 7a T2 

C8 ted 60 
C9] Ts, 

CUO den yh 

Ci Daght 

Cl voek20=GLGOxXC bt 038 
CS Date 

CVA Desi 

bios hs 

Cié6] u 

Biv he 


Ci82 A(2+JM)PZE+ CIM*CHF) 
Clea 7 


ea, te ’ y 
, : ) fe a x : F 
‘ ] + Tr 7 
; 7 ih o 1 
= an ; 
} in SON 
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; ray 7 
ae: wal 7 ny | 
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v 

Example 


To analyze a microstrip line with w/h=l and €&=9.9. 


The workspace SINGLE is loaded first: 


)LOAD SINGLE 


IMPD is executed next . 


Uae art ne Ue 
PEBECTIVE NIECECTRIG (CONS Pant 
&- OF 0198244 


CHARACTERISTIC IMPEDANCE 


91.2624 7005 


CHARGE DISTRIBUTIONS CH VS. ZE 


FIRST COLUMN? NORMALIZED DISTANCE 


SECOND COLUMN: CH 

OG PLS P11 S290] 5297746173934 
USS IVE LIS fsa oe ool 7d 
0.55358791487 S.854538571i 
0.7190733336 3-767108159 
0.8409902074 3.489595501 
OnFel6n07 2529) 3.6329 56052 
0.9682372823 3.598098021 
0.9905120983 5.580976234 
O. 77SSI9POSG? Ss a74G182535 
0.99994175 wei75g70186 


15s 


— 


eS 


‘\w thin endl ebtitedss fa) 8 cabana: * : fy 
| heed beat cana sa nay " 
. gxon petyoexe oh aay 


% 
4 


. 


ay ee 
inate be aaihaes 


p52 


FOTENTIAL DISTRIBUTION ON THE TAL 
FIRST COLUMN: NORMALIZED DISTANCE 
SECONT COLUMN? FOTENTIAL 


°26356032 0.641102785 
2+413403059 0.2778148304 
46596425771 0-.08741715987 
8.085810005 0-02427627784 


13.64080084 0-00750436892 


Appendix B 


Computer Programs for Coupled Microstrip Lines 

The workspace is called COUPLED. The parameters for 
numerical integration are the same as in SINGLE (App. A) so 
they are copied into COUPLED. Besides the symbols defined 
in App. A, there are some symbols used in COUPLED which are 


defined below: 


El through E8: AE) through [Eg]. Eqs 2-56 and 57 
FFl: [F,]; Eq. 2-58 

Gl through G4: [G,] through [Gy]; Eqs 2-50, 51 
Hl: ]H,] Eq. 2-58 

SIGs 6 Eq. 2-55b 

ZD: 1<@:<@' 

LE -1< G3< 1 

ZPi -< G5 <-] 

242 -OX< OKS 

cp he cee< 

ZPE SG je TERIGS el 

ZPF: ©), -0< Oi <-1 


ZP: a , aK or< So 
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Program PNT 


This program evaluates the coefficients and vector 


sizes and forms matrices which are necessary for the 


numerical integrations. 


v FNT 
Ci3 CF ¢6+2x#04 
C2] JDeEpnm 
CS JEe 9ZE 
C4] JF eeZF 
od JT¢ JD+JE+ JF 
C46] WWOe (ITs JD) PWM 
C7] WWEe CIT, JE) PW 
C82 WWFEC ITs JF) Pwr 


v 


Program DATA 


This program evaluates ¢; and q: for a given 6. 


v DNATA SIG 
Biel ZN€O0.5xX (SIG+1)+¢(SIG-1)xXNM 
C2) ZPOe C2xSIG)-Z0 
C34 ZPE¢(2xSIG)-ZE 
C4] ZPFe(2xSIG)-ZF 
Bi ZPeZPIts ZPEs ZPF 
C46] ZZ€Z0> ZEs ZF 


4 


Program INT1 
This program performs the integration of Eq. 2-53. 


v RREINTI 
Cid ALGTAICJI-THITtta 
ea AZCTH2ZCJIJ-THICII 
C33 R1ECCIAZ) &AQV=CTAL)*AL 
C4] BlLeCOL) —~fHLCIJISTHICT I 
C33 B2¢(01)-TH2CJI+THICTII 
C6 R2¢( CLBIO ABI) = GUR? ) ARS 
C7) C1¢(01)+THICJI+THICII 
C8) C2¢(01)+TH2CJI+THICII 
C9 RS¢CCI1C2)4C2)=+C ICL) aCl 
CiO0J R4e¢2x@R1IXR2xR3 
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Ci12 S$1i¢0.SxTHCJsI-THICII 

C122 S$2¢0.5xTHEJ#I]+THICII 

C133 S$3¢¢(10S1)+Si 

C141] S$4¢@(020S2)+((00.5)-S2)x(00.5)+S82 
CiSJ RSe€DTHEITX+/WMxS4 

Ci6] RR¢R4+RS-OTHC JIXCF 


v 


Program INT2 
Consider Eq. 2-50b 
Geaj =f. In(e\-6 /fr-g* de. (2-50b) 
Applying the Peed eeoenation 
¢=sing (2-32) 
to Eq. 2-50b yields 
. Ge, y=soq2in( ¢;' -sin@) *d@. (peat) 
Applying the identity 
(d/d@) [ (47/2) In( ¢, -sin@) ] = 
In( ¢,' —sin@) + 64/2) /( ¢; -siné). _ (B-2) 
to Eq. 2-50b results in 
Gg, y=2ln( |¢;' —sind,|*2-/*/|¢;' —sind ,|°1 7”) 
+2 f952[ (@-n/2) cosd/(¢;' —sin@] dé, ( B- 3) 
where the minus sign in Eq. B-2 is used. This choice of 
sign is used to remove the singularity occuring at ¢=1;, 
Considering that | 
G'; =2(0/n)- 4; (2-48) 
and the fact that ¢,; lies on the MIC for ¢& =l it is 
concluded that ¢;>1. 
For coupled lines with €,#1 an integral similar to 


that of Eq. B-1l occurs in which C. is replaced by ¢, and 


do?-S) ph Sp 2M poet sy em ie 
| ici carsctmngss uty 

ubtely dde-S pa a 
rt A (Ree Pye Ay a i a eee e 
cPhtqebs adit a ia A. 

| = { (tate idm Sa) CODNY, 

is) {Gate NAO > (hate OF i. 06 io " 
| ak situvean JOR Sep a 


(emit) Sb [Rates ! 
te soter= eta Seu &. 


isp Je gaenoOd.) ae aa 


rina 


(e-8) | 


can be less than or equal to minus one. Specifically, 
consider Hey ij ot Eqs. 2-57d, 58a, and 58b 
yay =e ims Ge (B-H) 
Applying the transformation of Eq. 2-32 followed by the 
identity of Eq. B-2 yields 
Hy, y>-2nin[ | sind gl "2*"/2/) ¢,-s ind [91477] 
—2nfo3?( 8+ /2) cos@/( ¢,-sind) dé. (B-5) 
where the plus sign of Eq. B-2 is used. 
Program INT2 evaluates both Eqs. B-3 and B-5. 


AZ can be ¢; of Eq. B-3 in which the minus sign of Eq. 
B-3 is relevant, or AZ can be Gs3/0t Eq. B=5 
in which case the positive sign applies. 


Correspondingly SGN can be -1 or +l. 


vy RS¢SGN INT2 AZ 
C1] SFe¢OSGNXO.S 
C2] R2¢(AZCII-10TH2C JT) eTHACII+SF 
C3] R1¢(AZCIIJ-1OTHIC JI) A THICII+SF 
C4) R3¢2xXx@IR2+R1 
CS R4e¢( (THEI I+SP)x20THC IJ¢$])=+AZCII-1OTHC Is] 
C6) RS¢OQTHEJTX+/WMXR4 
C7] RS¢ERS+KS 


v 


Program FIRST 
This program forms the matrix [G,] (Eqs. 2-51 and 58b) 


vy FIRST | 
Ci2 G1e(JE+JE) 90 
C2) et 


tow Cee er 

C4) L3sGilCilsJICINT1 
CSI Je J+1 

C6] 9CJZJE)/LS 

C7J ets 

C83 9CITSJE)/L2 


v 
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Program SCND 


This program forms the matricies [F,]» [G,], [H,], and 


\Eg\ (Basee235 lye 5G e etd. 58))c0 


v RACV SCND AZ 


Cid RAG (2kV) 20 
C2) I¢ei 
C3] L43Je1 
C4) LSSRACI#JI€VCZI INT2 AZ 
CSI JeJ+1 
C6] 9(J£VC2I)/L5 
C7] IeI+1 
C8) 9(T#£V01]3)/L4 
vw 
Program THRD 
This program assembles [G,] through |G]; {U} of Eq. 
2-51, [E,| through [Eg] and {F} of Eq. 2-56. The input 


to the program is VV, a vector which in this program contains 


Ey 


ie) 
hn 
C34] 
C4] 
C33 
C6] 
| 
C8] 
age | 
C1i0J 
Cit 
C12) 
C133 
C14] 
£15) 
C16] 
£17) 
£18] 
C19] 
C20] 
C21] 


ape eng ON. ro es 


vy THRE VV 

ER¢VVEL1I 

ET¢VVUC2) 

SIGeVVC3I 

DATA SIG 

ChLILEC2-ET) x2 

EL1¢((CER-1)x(SIG-1)=ET) xWWDO=+CF1+(ZZ°,-ZI) «2 
E2¢( CER-1)x(SIG-1)=ET) xWWO+CF1i+¢(ZFe,-ZIi) 42 
ES¢€(2x(CER-1)+ET) xWWF+CF14+(ZZ°,-ZF) 42 
F4¢(2x(ER-1)+ET) xXWWF+CF1i+(ZFe,-ZF)x2 
ES¢0O0.SXxWWEX@CF1+(ZZ°.-ZE) x2 
ES6¢O0O,SxWWEX@CF1+(ZFo.,-ZE)«2 
FFL¢(C IDs JEs 51) SEND ZD 

H1¢(JF>JE*1) SCNDI ZF 

E7¢(FF1+sC13 Gl)»f£1I3 Hi 

ESECIJTsUEs Ty SCN ZF 

Ut JEP 04 

G26C(JEs JE» 1) SCND ZPE 
G3¢-O0.S5xXxWWxeCFi+(ZEo.-ZE)4«2 
G4¢-00,5xXWWx@eCF1+(ZFE°,-ZE)*2 
GGE¢Gi+G2+G3+64 

GGO¢G1+G3-G2+G64 


v 


‘PHT LST vty 1) sdsoteram bdo enuot 
+ (8a ~ a 


te 10) ot Pmeoe ep ante monger. a 
(qn ad? «es so TS] te Fey duo 
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Program FNL 


The final assembly and rearrangement of the equations 


are performed by this program. 


EAL 


(i.e., to rearrange Eq. .2-56 into Eqs. 2-59). 


v FNL 
File” 300.5xXETX1-22 
F2€° 300.5XxXETX1+2Z 
F3¢~300.5XETX1-ZP 
F4e"300.5xXETX1+ZF 
FE@(C1-ER)XFItPetrs+F 4 
FO¢(1-ER) XFItF2-FS+F4 
CF2¢02X1+ER 
VieEJIit i JE 
PEEVIJCFECVII-Cr2 
-FOCVIIJ¢FOCVII-CF2 
GES CE1+E2) +» (E7t+ES-ES+E46) »ESt+E4 
GO¢(E1-E2)>(E7-E8+ES-E6)»E3-E4 
Fer 
L6é:;GECI¢IJ¢GECIs#II+CF2 


£157 GOCIsIJ¢GOCI+sII+CF2 
ELiG) Tel+1 
C172 2cCIHJD)/L6G 
Grew? - fel +e 
£19) L7Z3GECIsII¢GECI¢II+CF2 
£20] GOCIsIIJ¢GOCIs1I13+CF2 
Gers Ferri 
B22 y) (9 CI S2J5T)/L7 
vv 
Program FNL2 


are solved and the various coupled mic line parameters are 


The final equations for both cases 


€p=l and € 41 


evaluated by this program. 


Crd 
C2] 
C3) 
C4) 
CSI 


vy FNL2 
ABCE+¢FEBGE 

AE¢ JII*ABCE 

BE¢ JDY(-JF)+ABCE 
CE¢(-JF)fABCE 
ABRCO+FOEGO 
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C6] 

C73 

C8) 

C9IJ 

£10] 
BLRID 
C12) 
£13] 
C14] 
Ea 
C16] 
ey By a 
£18) 
Oi7 J 
C20] 
C21) 
C22] 
C23] 
C24) 
C25] 
C26] 
C27] 
C281 
C29) 
C30] 
C31. 
E32) 
£331 
C34] 
C35) 
C36] 
C37) 
C38) 
a 
C40] 
C414 
C42) 
C43] 
C44] 
C451] 
C46] 
C47] 
C48] 
C49) 


ADe JIIFABCO 

RO¢ JI (- JF) VABRCO 
CO¢+(-JF)tABRCO 
CE1¢+/WJIXBE 
CO1¢+/WJIXRO 
CEQ0¢+/WJIXUBGGE 
CO0¢+/WIxXUEGGO 
T12 

ae 

T1 


240+(C ICEOXCE1) 40.5 


T12 
T12 
T2 


240+¢€1COO0xXC01)40.5 


T12 

fie 

pea 
ICE1=CEO 
TL 

Ti2 

TSP 
{C01=+COO 
T12 

Ti2 

1S § 

T4 

TS 

TS 
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T10 


R(S*JIF)9ZF+CEsCO 


v 


159 


Lae 
wei 
i 
‘ 


a = 
i, 
/ J 
iy a 
a : 
s15 
, 
i - 
a4 
ne 
~ q — 
‘ 
~e 


s 
Co 
<¥ 
© 
/ 
7m , 


i” 
ui 7 , 
: iy 
i F) 
7 

( 

‘ 


160 


IMPD 


This is the main program. The input to this program 


is VV. In this program VV is a vector with the following 


elements: 
VV[11 =ér, 
Ww (2|=w/h (Fig 2-10a) and 
VV[3]=s/h (Fig 2-10a). 
VV is modified by IMPD to contain &,n , and S' needed as 


the input to THRD. 


v IMFD VV 
ie VVC3IE1L+VVUCSISVVUL27 
f23 VUC2I¢VVLC29<+2 
C3] THRD VV 
C4] FNL 
CSI PNEZ 


v 


Program TTL 


This program produces the titles and comments necessary 


for the output in FNLé2. 


v TTL 
ee Be T1l¢*EVEN-MOLIE CHARACTERISTIC IMPEDANCE! 
C2] T2¢'OND-MODE CHARACTERISTIC IMFEIIANCE' 
C33 T2P¢‘EVEN-MODE EFFECTIVE [DIELECTRIC CONSTANT! 
C4] TSP¢*ODUD-MODE EFFECTIVE DIELECTRIC CONSTANT' 
CSI T3¢*CHARGE NISTRIBUTION? ' 
C6) T4¢'FIRST COLUMN: NORMALIZED DISTANCE ON THE MIC’ 
C73 T5¢'SECOND COLUMN? EVEN-MOIE CH' 
C8) T6€' THIRD COLUMN? ODD-MODE CH' 
Eos T7€*FOTENTIAL DISTRIRUTION ON DAI RETWEEN MICS' 
[107 T8e¢'FIRST COLUMN: NORMALIZED DISTANCE' 
Ci1d T9¢*SECOND COLUMN: EVEN-MODE FFROTENTIAL' 
C123) T10¢'THIRE COLUMN: ODD MONE FOTENTIAL' 
C13] Tile POTENTIAL DISTRIBUTION ON DAI OUTSIDE MICS! 
Cina? ieee 
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Example 
To analyze coupled mic lines with 
Css | 
w/h <=1.0 , and 
s/h =1.0 


the workspace COUPLED is loaded first. The response to 


IMF 9 a1 


is 


EVEN-MODE CHARACTERISTIC IMFPETIANCE 
7% /0066932 


ODU-MODE CHARACTERISTIC IMPEDANCE 
43.2229°4308 


EVEN-MODE EFFECTIVE DIELECTRIC CONSTANT 
7 ¢856033191 


ORUCMOUE EFFECTIVE NIERECTRIC CONSTANT 


Je 22740583 


CHARGE TISTRIEUTION: 

FIRST COLUMN: NORMALIZED DISTANCE ON THE MIC 
SECOND COLUMN: EVEN-MODIE CH 

THIRE COLUMN: OFT-MODE CH 


“Oo9PS994AL7ZS 6¢ 794555263 7 eo PE4AGI77IO4 
"O68 9SS7I0GG?  SePaslGa7e4 6.8120589 

“O04 9790S51209SS S.6235646SS Oe S591 OL 607 
$0. 9682372823 3S. S706509%S 6-003415031 
(ORF eUGar ages, Ose SISO 4e2 2° 799411158 
TO. 8907702074 32158074753 3+6703186S8 


32984286331 
7 e OD4407753 


We71I9OFS3S55 
TOvooses 7 LAS s 


7° 634514068 


UPOOT OAD? LG 


& zeny Oro? 
Re nnd eee 


J eVeLe lee Se 


we 7IO/7127785 


TO.1199182902 4.959754888 3e«886318488 
O.119°9182902 4.848171275 5 -9042027748 
0.4502192597 4.4688702156 66224333422 
0.5538791487? 4.506285863 6+435227221 
0.71907 33836 4.342177833 6<672403318 
O.SMOSSO2074 4.239552 3542 6 - 928322168 
O.PEP6992529 A. 229550085 7198439857 
0.96S82372823 4.3528538908 7¢491717851 
O.8905120983 4.549176108 7 840944982 
O.F9ESS9O0SS9 4.F009S1007 8. 289778912 
0.99°94175 Jee? 7188024 92 SS5116178 
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FOTENTIAL DISTRIBUTION ON DAI OUTSINE MICS 
FIRST COLUMN: NORMALIZED DISTANCE 

SECONI COLUMN: EVEN-MODE FFOTENTIAL 

THIRD. COLUMN: OLD MODE FOTENTIAL 


~1.26356032 0-6779665161 0.6687968417 
“24413403059 0.3032838761 O.2759921618 
“4.596425771 0.101707095 0.07841470369 
~8.085810006 0.03033464862 0.01752825718 
~13-.64080084 0.008757500631 0.003875702662 


POTENTIAL NISTRIBUTION ON DAI BETWEEN MICS 
FIRST COLUMN? NORMALIZED DISTANCE 

SECOND COLUMN: EVEN-MOTIE FFOTENTIAL 

THIRD COLUMN: O00 MODE FOTENTIAL 


1.0130446736 0.99558 S6221 0-9516569076 
1.057458317 0.906074752 0.7938957452 
te lLSGO2L9S216 0.5321237625 0.6461210026 
1 2835302505 OV 7G So 7 SS 0.509282825 

16425562831 0.7208040195 0.3846162811 


| -taaw bay, 4 
1.5744371659 
=< bay J . : 

Le PIKES TOO? 
1.839704784 


1.932531 683 
1.9669 53264 


0.48453872529 
0+ 553052949 
O.G54 3171934 
0 .6501500673 
0. 64929249711 


0+-27345462835 
0.177495482 
OCP RLS SI 7E9? 
0.04152114196 
0.008021078655 
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Appendix C 


Workspace THAT 


This workspace contains several functions for evalua- 
ting and solving the final equations of Chapay le. The 


notation used is as follows: 


B, BB: B,Eq. 4-234, 
Cr C,Eq. 4-23e, 
DEL 6,Eq. 4-222, 
EF: Ce, 
: an 
ET: y,Eq. 4-22c, 
: (fh) normalized frequency CHz-mm 
Gliltls Jdtgeeane Gig overs Eqs. 4-33 
Ls TE Edi t— OL 
Ms M,Eq. 4-32a and 32c 
PH=tan"(C/L) C,and L are defined in Eqs. 4-23e and 31 
RO: h ,Eq. 4-25, 
SG: On,Eq. 4-25, 
WI [ I]: W; , Fig. 4-3, 
>a bel GE é,,Fig. 4-3, 
meliys | 0 Wig $13, Fig. 4-3, 
XU [I]: ae anges EET 
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C4] 


FO computes ei; and Su given the values of ee and Wy - 


vy FO WI 
Nepwl 
XUeENPO 
J¢O 


Pi beJeg+i 


XUCJIJe+/ ItWI 
ACISNO/LIL 
XL¢0O»% ~“1¥XU 


v 


Dl computes the matricies with elements 


(Eq-fd. (éy€d. (Eqtfst2n), and ( Eyt€st27) 


used as the arguments of Eqs. 4-33a, c, and d. 


v Dl 
XIT1¢XI°.-XL 
XI2¢XIe.-XU 
Z1ieXIe.t+XL+2xET 
Z2¢X1°.+XUt+2XET 
MePxXxl 


v 


D2 computes the factors 


exp(—€aL), exp(—-fyL), exp(-€,L), 


which appear when integrating Eqs. 4-33. 


Eid 
C2] 
C3] 
C4] 
ES2 
C6] 
E74 
C8] 
Cot 


vo BAIL 
EI¢(M+M) 90 

I¢0 
LiO?I¢I+1 
EICI¢IJ¢a-LxXICII 
2(I<M)/L10 
EX1L1¢(M»M) 9 *-LXXL 
EX2¢ (MoM) PA-LXXU 
LX2¢ (MoM) PLXXU 
LX1¢(MsM) PL XXL 


Vv 
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D3 computes the matrices SNU and SNL such that 


correct signs can be allotted to the different terms for 


evaluating the integral in Eq. 4-33c, for the given conditions 


Cid 
C2) 


oy 


forms matricies Gl, G2, and G5 with elements G 


Geiay 


Cid 
bed 
C3] 
C4] 
C33 
C6] 
C7J 
C83 
Cea 
C10) 
Ciij 
Ci2] 
Ci3si 
Ci4] 
Lisa 
Cié) 
C17) 
Cis 
ivi 
C20] 


v DS 
SNUC2xX7~0.5+XL°.<xI 
SNL¢2x~0.5+XUe.>XI 


v 


Using the matrices computed in Dl, D2, and D3, FL 


a3)» 


» and Gs ij respectively. 


ey OR ax ae 

Se(La22)4+Ca2 

SQ¢S40.5 

CLI¢CHL 

FHe~ 3oOCLI 

LIQ¢e (FH+CLI)=L 
E2¢EX2=+2xSQ 

EL1¢EX1+2xSQ 

RiI¢PH-CxxI1 

R2eFH-CxXI2 
R3¢-SNUXE2x(LIQ+XTI2)x(10R2)+R2 
R4¢-SNLXE1L1X(LIQ+XI1)x(C LORI )=+Ri 
G2eRZS+R4+EI=FS 

RS¢PH+CxZ1 

R6¢FPH+CXZ2 
R7Z¢ELx(LIQ+Z1)xCLORS)=R5 
R8¢E2x (LIQ+Z2)x(10R6)=K6 
Git¢R7-R8 

R9¢E10C1s JX2OPHtCXET+XL 
RLIO¢E20L13 7 x20FH+CXET+XU 
GS¢R9Y-R10 


v 
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tt . 
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FS 


matrix equation and evaluates the determinant of [G] (Eq. 


4-35) ; 


F2 finds the matrices G3 and G4 (Eqs. 4~33a and b). 


vy ESFe C 

G4¢G63¢0 

N¢O 
LiIsNeNti 

RS¢( (ON) 42)-Ca2 

ROERS*0.5 

DR€=+ROXRO+L 

GS¢G3+tDRx (a-LX1+ROXZ1) -4-LX2+ROXZ2 
G4ieEX1x¢a-ROXIXI1) 
641¢G41=ROxXL-SNLXRO 

G42¢EX2x Ca-ROXIXI2)=ROXL+SNUXRO 
G4I¢2xEIe(La2)-RS 
G4¢64+641-G642+64I 

9OCN<MAX)/L1I 


v 


similarly, F3 finds’ G6 and G7 (Eqs. 4-33f and 2); 


vy L FS B 
66¢G67¢0 
Nev d 


L23NeNt+i 


SGS¢€(BA2)+(ONFDEL) &2 

SG¢SGSr0.5 

BRSG¢R+SG 

DR¢ RSG+SG+tL 

G6e¢G6+t0RX (C4-LX1+SG6XZ1)-*-LX2+SG6XZ2 
G71i¢+BSGXEXix¢a-SGXIXI1)=+L-SNLXSG 
G72 BSGXEX2x (a-SGXIXI2)+L+SNUXSG 
G7I¢2xEI=(L4’2)-S6S 
G7¢G67+G671-G72+G67I 

9CNMAX)/L2 


v 


Having the necessary input, F5 assembles the final 


DET is a function for evaluating the determinant of 
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a matrix. 


Cid 
C2] 
C3] 
C4] 
Cod 
C6] 
C7J 
C8I 


Fixing the frequency, F5 should be executed few times around 


the expected € , until the point is found for which the 


VmrE SS ER 
CeCixFx(ER-EF) 40.5 
RRECLixFXx(EF-1)40.5 

Leet Cc 

ar 2.0 

L F3 BE 

G¢G1+G2-G63+G64 
HeGSsCLICERXDELXBBXG)-G6+G7 
NET 14 14 #H 


v 


determinant of [G] vanishes. 
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Appendix D 


Workspace RTI 
The following notation is used in this workspace: 
EPS: the precision to which the roots of Eq. 5-6c 
should be found, 
BRel + 46, 
Fs fh in GHZ-mn, 
NX number of positive roots of Eq. 5-16b to be found, 


RO: An 


fe) 


This function evaluates Eq. 5-16b. 


ev RReED Ri 
Pe) RS¢Ris2 
Heil MUeC1LAS-RS)*0.5 
fa 4CASERS)/LI 
C43 FiemMUXSONELxXmMuU 
Baal F2¢s6onEL xMU 
oud = 3 a 
E7ie? CVS Te-MUXTORELXAU 
CSI F2¢20DEL xMuU 
C9J L2tRR¢ECERXFIXx20R1)-RIxF2x10rR1 


v 


ROOT 
Given a certain frequency, ROOT finds f,» the TITSseeNnn 


positive roots of Eq. 5-6c. 
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ciel 
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ee 
C30] 
si. 
Goce 
C33] 
C347 
Bea | 
GSa) 
a7 


vy F ROOY NX 
AS¢(ER-1)x(CCLixF dae 
RO¢NXPO 
N1¢N¢O 
LOS NEN+) 
FOCNSNX) 7/0 
Lisl¢o 
AtO0.SXN+N1-1 
KteOO.,SXN+N1 
L23MeO 
VeAt+TXE-A 
LSimMens1 
Lett 
ado Rea ye a Se 
RSi¢D VCMI 
RSet le VL 
SNERSLIXRS2 
9(SNO0)/L4 
9(SN=0)/L5 
>OMS1.0) ALS1 
ay yt 
LSI 3NIeCNI +1 
sc a 
L4:tAtVEMI 
Ke VCM+tI I 
4CCR-Ad SERFS) /Le 
ROCNI¢VE MI 
ILO 
LSt3CRSL=0)/LS 
ROCNI¢VCM+iI 
alo 
L&é7;ROCNI¢VCNI 
ILO 
L73'NO CONVERGENCE ' 
PN CIS" 
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‘ITERATION NO. IS' 
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Workspace NIN 


The notation used in this workspace is as follows: 


As is A, Eqns o-7C+ 
Cs is Cy. Ede 5-765... 
EF: 1s &) 
: is &) 
ET : is, Eq. 4-22c, 
M1: is M, Eq. 5-10a, 
M3: is the number of terms retained in Eq. 5-10a, 
s sR, fd. 5-15; 
NDB + is %, Eq. 5-9c, 
NMB : Steves Eq. 5-9c; 
RO: Pm? 
TH: | Ges Fige 5-35 
THL : Oy: 9 Fig. 5-3, 
THU : 6,35 Fig. 5-3. 
DATA 


Finds Oj.) @., and 6,, (Fig: 5-3). 


v DATA 
es Ee MVe pl 
C2] THUtO0.5x (MR) =MR 
C3 THL¢O+ ~14¥THU 
C4) SLELOTHL 
C5] SU¢1LOTHU 
CéJ SM¢O.SXSL+SU 
C74 tice Loony 
C8] MER¢O.SXTH-THL 
aR FREO.SXTH+THL 
C101 MIR+¢0O.SxXTHU-TH 
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Computes f, of Eq. 5-38b. f) of Eq. 5-20a is included. 
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Computes f3 of Eq. 5-20c. 
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Computes the coefficients of Eq. 5-10a using the roots 


of Eq. 5-6c; 2 should be copied from RT first. 
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FH is executed around the expected values of € for a 
eertain fh(tor which fs found in RT) untilvan. €, is 


found for which the determinant of JH] vanishes (Eq. 5-16). 
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» FH finds the determinant of [H] using DET. 
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